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FOREWORD 


This  report  describes  the  work  performed  by  the  Cornell  Aeronautical  Labora¬ 
tory  under  Contract  Nonr-3659(00)  (FBM),  sponsored  by  the  Office  of  Naval 
Research  of  the  Department  of  the  Navy*  The  time-period  covered  is  from 
16  October  1961  to  25  March  1962.  The  program  is  being  conducted  under 
the  general  direction  of  Cdr.  F.  R.  Haselton,  Code  466,  and  R.  Cooper, 

Code  438,  of  the  Office  of  Naval  Research. 

The  authors  wish  to  express  their  gratitude  to  Mr.  D.  C.  Clark  of  the  Vehicle 
Dynamics  Department  of  CAL,  for  his  contributions  to  the  stability  and  con¬ 
trol  studies,  and  to  Messrs.  Q.  Tufts,  E.  Sullivan,  and  F.  DuWaldt,  of  the 
Applied  Mechanics  Department  os  CAL,  for  performing  the  theoretical 
hydrodynamics  studies  and  for  valuable  contributions  to  the  entire  program. 
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ABSTRACT 


Hydrodynamic  and  stability  and  control  characteristics  are  investigated  for 
a  novel  submarine  configuration  employing  tandem,  large  hub-to-diameter 
ratio  propellers  whose  blades  are  pitched  collectively  and  cyclicly.  Stability 
coefficients  for  the  hull  and  the  propeller  are  derived.  Expressions  for  the 
control  forces  and  moments  generated  by  the  propellers  are  developed  and 
the  perturbation  equations  of  motion  are  used  to  analyze  the  controlled 
dynamics  of  the  submarine.  An  analysis  of  the  motions  in  the  ve-tical 
plane  indicates  that  the  Submarine  can  be  stabilized  at  high  speeds  by  auto¬ 
matic  control  of  collective  and  cyclic  pitch  and  that  reasonable  maneuver¬ 
ability  can  be  obtained. 
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I 

INTRODUCTION 


In  October  1961  the  Cornell  Aeronautical  Laboratory,  Inc.  undertook  to  inves¬ 
tigate  the  hydrodynamic  and  stability  and  control  aspects  of  a  novel  submarine 
configuration  utilizing  variable-pitch  large  hub-to-diameter- ratio  propellers. 
The  configuration,  invented  by  Cdr.  F.  R.  Haseiton  of  the  Office  of  Naval 
Research,  employs  two  propellers  mounted  circumferentially  (forward  and 
aft)  on  a  neutrally  buoyant  body  of  revolution,  to  produce  any  combination  of 
forces  and  moments.  By  means  of  this  arrangement,  which  has  been  called 
a  Tandem  Propeller  Submarine  (TPS),  it  is  possible  to  produce  control  forces 
in  all  three  degrees  of  freedom,  or  control  moments  in  all  three  degrees  of 
freedom,  as  well  as  combinations  of  forces  and  moments. 

The  work  undertaken  by  the  Cornell  Aeronautical  Laboratory  has  as  its  basic 
objectives  (a)  the  theoretical  determination  of  the  hydrodynamic  characteris¬ 
tics  of  the  tandem  propeller  configuration,  (b)  an  investigation  of  the  trim 
and  stability  and  control  characteristics  of  the  controlled  submarine,  and  (c) 
a  comparison  between  the  TPS  and  a  conventional  submarine  with  respect  to 
stability  and  control  characteristics  and  handling  qualities. 
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The  work  accomplished  to  date  include*  the  following: 

(1)  The  hydrodynamics  have  been  defined  in  terms  of  hull  and 
propeller  stability  derivatives  for  the  high-speed  case  (about 
24  knots).  Work  on  the  hovering  caso  is  in  progress. 

(2)  The  equations  of  motion  for  the  high-speed  case  have  been 
written  and  numerical  values  for  all  of  the  physical  and  hydro- 
dynamic  properties  have  been  computed  for  the  postulated 
submarine  configuration. 

(3)  The  high-speed  trim  characteristics  have  been  investigated 
for  operation  with  one  propeller  and  operation  with  two 
propellers. 

(4)  The  propulsion  and  control  coefficients  defining  the  generation 
of  forces  and  moments  as  a  function  of  collective  and  cyclic 
pitch  control  inputs  have  been  determined  for  the  high-speed 
case. 

(5)  An  analog -computer  study  of  the  pitch  plane  dynamics  of  the 
controlled  (i.  e.  stability  augmented)  submarine  has  been 
performed  for  the  high-speed  case.  The  effects  of  pitch 
maneuvers  on  roll  angle  and  forward  speed,  as  well  as  the 
effects  of  certain  nonlinear  control  coupling  forces  and  moments 
on  pitch  dynamics  have  been  determined. 

(6)  Tentative  control  system  gains  for  pitch  attitude  and  depth¬ 
keeping  loops  have  been  selected.  Two  types  of  stabilization 
have  been  investigated  and  their  differences  noted  for  the 
pitch  maneuvering  mode  of  operation. 

(7)  A  limited  analytical  study  of  the  yaw  plane  dynamics  has  been 
accomplished. 
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STATUS  OF  PROGRAM  It  SUMMARY  OF  CONCLUSIONS 


2.  I  STATUS  OF  PROGRAM 

(1)  High-speed  hydrodynamics  are  defined  for  the  TPS 
configuration,  permitting  a  high-speed  stability  and 
control  investigation  to  ascertain  the  central  issue  of 
feasibility,  namely:  is  sufficient  cyclic  pitch  available 
to  stabilise  and  provide  maneuver  control? 

(2)  Low-speed  (or  hovering)  hydrodynamics  are  tentatively 
defined  for  both  the  propeller  and  hull  contributions. 
Calculations  must  be  performed  to  assess- 

(a)  The  relative  contributions  of  the  tandem  propellers 
and  the  hull  to  the  total  hydrodynamic  forces  and 
moments  acting  on  the  TPS. 

(b)  The  degree  to  which  the  propeller  forces  are 
simultaneously  influenced  by  combined  blade* 
pitch  changes  and  changes  in  blade  flight  path-angle 
due  to  the  motions  of  the  submarine. 
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(3)  Equilibrium  trim  studies  have  been  made  to  ascertain  the 
ability  to  trim  the  TPS  in  steady  longitudinal  flight  when 
one  propeller  only  la  used  for  propulsion. 

2.2  SUMMARY  OF  CONCLUSIONS 

(1)  Analysis  of  the  stability  coefficients  for  the  high-speed 
case  shows  that,  in  comparison  with  conventional  sub¬ 
marines,  the  gains  in  control  effectiveness  achieved  at 
low  and  zero  speeds  {plus  the  gain  in  control  flexibility) 
are  obtained  at  a  sacrifice  of  high-speed  control  effective¬ 
ness. 

(2)  The  TPS  is  dynamically  unstable  in  pitch,  unless  control 
forces  and  moments  are  applied  to  modify  a  id  eliminate 
this  instability.  The  instability  arises  primarily  from  the 
large  unstable  pitching  moment  due  to  angle-of-attack  (//„  ) 
relative  to  the  levels  of  damping  in  pitch 
by  the  TPS. 

(3)  The  divergence  in  pitch  motion,  due  to  instability,  is 
sufficiently  severe  to  require  automatic  stabilization  in 
contrast  to  manual  stabilization. 

(4)  Although  other  types  of  stabilizing  feedbacks  were  investigated 
for  use  in  an  automatic  pitch-control  system,  simple  pitch- 
rate  and  pitch-angle  feedbacks  were  found  to  be  effective. 
Well-damped  pitch-angle  responses  to  pitching  moment 
control  inputs  can  be  obtained  in  about  60  seconds.  Corres¬ 
ponding  steady-state  dive  or  climb  rates  of  about  6  ft/sec.  are 
easily  achieved.  These  figures' do  not  represent  maximum 
achievable  pr rform.i nee. 
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(5)  An  analog  compute) r  investigation  shows  that  the  lnvni  of 
the  command  control  inputs  must  be  limited  to  keep  the 
propeller  blade  anglus-of-attack  within  their  stall  limits 
when  cyclic  and  collective  pitch  are  used  for  automatic 
stabilisation.  Computed  feedback* loop  gains  that  originally 
appeared  to  be  so  large  that  continual  saturation  of  the  eye  lie  - 
pitch  control  would  result,  were  found  to  be  acceptable,  pro¬ 
vided  limits  were  placed  on  maneuver  demands, 

(6)  Although  an  investigation  of  the  yaw-plane  dynamics  has 
not  been  completed  it  is  probable  that  conclusions,  similar 
to  those  given  above  for  pitch,  will  eventually  be  reached. 

(7)  It  has  been  shown  that  within  certain  speed  limits  it  is 
possible  to  trim  the  TPS  in  high-speed,  straight  and  level 
flight,  with  one  propeller  fixed  and  one  operating.  Speeds  of 
about  15  knots  can  be  achieved  with  a  single  propeller  operating 
at  50  rpm.  (An  investigation  of  the  problem  of  maneuvering  the 
TPS  with  one  propeller  operating  remains  to  be  completed.  It 
is  anticipated  that  this  control  mode  will  present  a  serious 
problem.  ) 

(8)  With  respect  to  the  overall  stability  and  control  prcblems,  the 
tentative  conclusion  is  reached  that  automatic  control  of  the 
TPS  submarine  is  feasible.  Comparison  with  a  conventional 
submarine  has  not  yet  been  made  nor  has  the  question  of 
handling  qualities  been  examined  fully. 
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SYMnOLS  AND  NOMENCLATURE 


3.1  NOMENCLATURE 


Forces  and  moments  acting  on  the  Tandem  Propeller  Submarine  (TPS)  are 
written  in  body  axes.  The/'X-y-^  body  axis  system  is  a  right-handed  ortho¬ 
gonal  triad  with  its  origin  at  the  submarine  center-of-g r  vity.  ft*  is  positive 
forward,  is  positive  to  starboard  and  y  is  positive  down.  Components  of 
total  linear  and  angular  velocity  along  and  around  ftt  -y- y  are^./l/",  MT  and 
p  t  ,  si*  respectively.  Generalized  forces  and  moments  are  ^  -  Y  -  ?  and 
K  -  M  -  N  along  and  around  /P  -  -  y  . 


Trim  conditions  will  be  denoted  by  a  zero  subscript  and  properties  related  to 
forward  and  aft  propellers  will  be  distinguished  by  subscripts  ~  and  . 
Perturbations  from  trim  m  velocity  are  denoted  by  barred  quantities,  thus; 


-  n/'  +  nr  ,  ,ur  -  A/r.  s.tr  (l> 

*  >  j 

Angular  velocity  perturbations  are  defined  by: 

(p  ‘  fa  +  $  ,  f  ••  fo  *■  *>  ,  si  ‘  SI*  ■*-  'P  O) 

"her*'  ^  ,  O  and  ^  .ire  the  linearized  (small -.ingle)  Euler  angles. 
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Stability  derivative#  wilt  be  identified  by  subscripts  which  indicato  the  perturbed 
variable.  For  example,  X#  is  an  /#  -force  derivative  or  coefficient  denoting 
a  change  in/#  -force  duo  to  the  -velocity  perturbation  .46  , 

Other  symbols,  assumptions,  etc,  will  be  presented  as  they  occur  in  the  text. 

In  addition,  Section  3.2  contains  a  list  of  all  of  the  symbols  used  in  this  report. 


3.2  SYMBOLS 

All  symbols  used  in  this  report  are  listed  below.  Where  used,  the  subscripts 
£  and  66  denote  forward  and  aft  conditions. 


Hull  Characteristics 


Units 


dL  - 

Max.  hull  diameter 

ft 

f .. 

Fineness  ratio  i 

dimensionless 

i  - 

Length  of  submarine 

ft 

A- 

Distance,  fore  point  to  c.  g. 

ft 

L- 

Distance,  prop,  plane  to  c.  g. 

{  carries  no  sign) 

ft 

/•  = 

Metacentric  height 

ft 

<4  = 

Submerged  displacement  (=  Weight) 

lbs 

Propeller  Characteristics 


R  , 
D  , 


Average  radius 
Average  diameter 


ft 

ft 


N  :  Number  of  blades 

-  Blade  area  *  tt 

-Q  -  Blade  angular  velocity  (  Cl  carries  no  sign)  rad/ see' 

TL  -  Blade  angular  velm  ity  (  Yl  carries  no  sign)  rev/jrtec 


| 
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Unit* 


Lift  coefficient 

Q.’ 

Profile  drag  coefficient 

//  - 

Induced  drag  coefficient 

T  : 

Flight  path  angle  of  the  blade  (see  Fig.  1-2) 

rad 

: 

Blade  angle -of -attack  (see  Fig.  *1*2) 

rad 

cf  - 

Total  instantaneous  blade  pitch  (see  Fig.  '1*2) 

rad 

<£: 

Collective  pitch  (see  Fig.  -1-2) 

rad 

4>cf, 

Change  in  collective  pitch  (see  Fig.  4-2) 

rad 

4  •• 

Blade  angle  of  attack  due  to  collective  pitch 
(  b)  :  <£>  *•  A'f  ~  ) 

rad 

cr  z 

Blade  a/umuth  angle  in  piano  of  prop 
(see  Fig.  1-2) 

rad 

6^  (7^  ; 

Sine  component  of  cyclic  pitch  (see  Fig.  4-2) 

rad 

c/^c^  ^  - 

Cosine  component  of  cyclic  pitch  (see  Fig.  4-2)  rad 

Kinematics 


U-- 

Total  velocity  of  c.g.  (  (J  ^  M."*"*-  sw^rstv'*' ) 

ft/  sec 

M-  - 

Ol  -component  of  velocity  (perturbation  =  AL  ) 

ft/  sec 

nr. 

-component  of  velocity  (perturbation  =  OS’) 

ft/  sec 

ACT  - 

-component  of  velocity  (perturbation  -  ) 

ft/  sec 

P* 

f/'  -component  of  total  angular  velocity 
(perturbation  --  <fi  ) 

rad/ sec 

If¬ 

V-  -component  of  total  angular  velocity 
^perturbation  &  ) 

t*ad/  sec 

ni  , 

-component  of  total  angular  velocity 
^perturbation  .  <p  ) 

rad/ sec 

V 

Blade  velocity  relative  to  water  (see  Fig.  1  2) 

ft/  sec 

Mass/lnortia  Properties 


Units 


KrKi-Kj* 

Coefficients  of  accession  to  mass  adong 

dimensionless 

/<*-*,-<• 

Coefficient  of  accession  to  inertia  around 

dimensionless 

m  = 

Mass  of  submarine 

#-sec  / ft 

m,= 

Virtual  mass  along =  /7£  (/  +  Kt) 

# -sec2/ ft 

Virtual  mass  along  ^  =  771  ( f  + 

#  -sec2/ft 

TTij  = 

Virtual  mass  along  y  =  (/  *•  Kj ) 

#  «see2/ft 

Submerged  moments  of  inertia  about 

„  ,  2 

#  -ft-sec 

L* 

Virtual  moment  of  inertia  about  -  T^( '  +K?) 

#  -ft-sec2 

4/  * 

Virtual  moment  of  inertia  about  <y*Iy(l  +  Kf) 

#-ft-see2 

«r- 

ii 

Virtual  moment  of  inertia  about  y  ■  Iy  (f  +  Kt) 

#-ft-sec2 

H  , 

Propeller  angular  momentum 

#-ft-aec 

e  = 

Water  density 

# -sec2/ ft4 

Forces  and  Moments 

^X.0  = 

Trim  thrust  available  from  propeller(s) 

if 

Aco» 

Trim  moment  at. the  propeller(s) 

ft-# 

Po  = 

Trim  propeller(s)  power 

hp 

'V 

Metacentric  pitching  moment  coefficient 
*  -”1 1 /» 

ft-#  /  rad 

Vl'cV 

Control  (propeller)  forces 

# 

KcMc"c- 

Control  (propeller)  moments 

If  ft 

9 
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.■ 

zx-sr-zz* 


ZK-Lm-ZN* 


Combined  propeller  and  hull  hydrodynamic 
and  hydrostatic  forces 

Combined  propeller  and  hull  hydrodynamic 
and  hydrostatic  momonts 


Unit 

I 

#-ft 


Hull  Derivatives 


X 

Aajl  = 

Dimensionle 

11 

Dimcnsionle 

Dimensionle 

•ct 

II 

Dimensionle 

Dimensionle 

If 

*  ^ 
fll 

Dimensionle 

II 

Dimensionle 

N* 

II 

Dimensionle 

VI 

y 

Dimensionlei 

Dimensionlei 

w'f  - 

Dimensionlei 

N!r  = 

Dimensionlei 

K- 

Dimensionlei 

N'f- 

Dimensionlei 

sfr  -force  derivative  due  to 

fj>  -force  derivative  due  to  /V" 

0 

-force  derivative  due  to  j*' 
-force  derivative  due  to  $ 
f  -force  derivative  due  toxj^ 

^  -force  derivative  due  to  3 
-moment  derivative  due  to  AT 
^  -moment  derivative  due  to 
-moment  derivative  due  to  ft 
y  -moment  derivative  due  to  MA 
-moment  derivative  due  to  b 
]r  -moment  derivative  due  to  AT 

m 

f  -moment  derivative  due  to  p 

9 

-moment  derivative  due  to  fi 


Propeller  Terminology 


vX 
$  ^ 
II 

Dimensional 

—  * 

&  -force  derivative  due  to 

#/ft  per  sec 

xp 
y  = 

Dimensional 

&  -force  derivative  due  to  $>p 

Hi  rad  per  sec 

YP  - 

•at 

Dimensional 

p  -force  derivative  due  to  Ar^ 

#/ft  per  sec 

*  Tho^  subscript  denotes  perturbations  with  respect  to  propeller  axe h , 
Propeller  axea  are  parallel  to  body  axes  but  with  origin  at  x  >  , 

forward  and  aft.  ”  ~ 
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y 

't  ■ 

Dimensional 

Unit 

f  -force  derivative  due  to  Qp  # / rad  per  sec 

Dimensional 

^  -force  derivative  due  to  #/ft  per  sec 

iX  . 

Dimensional 

-force  derivative  due  to  #  /  rad  per  sec 

*/■ 

Dimensional 

^  -moment  derivative  due  to  ^  # -ft/  rad  per  sec 

KL- 

Dimensional 

-moment  derivative  due  to -5^  # -ft/ ft  per  sec 

nf. 

Dimensional 

-moment  derivative  due  to  Op  # -ft/ rad  per  sec 

Dimensional 

^  -moment  derivative  due  to^  # -ft/ft  per  sec 

K* 

Dimensional 

^  -moment  derivative  due  to  M»y  # -ft/ft  per  sec 

NX  -. 

Dimensional 

^ -moment  derivative  due  to  fp  # -ft/  radper  sec 

F*-Ft-Fy  - 

components  of  lift  and  drag,  # 

prop,  axes 

Ft  -- 

Tangential  force  at  average  radius  R  # 

components  of  moment,  prop,  axes  #-ft 

ft  ■ 

Trim  thrust, 

moment  and  horsepower 

i.-£>  = 

Lift  and  drag  forces  ^ 

n 

N> 

<J 

X 

Dimensional 
due  to  /a  6 

tf-  -force  propeller  coefficient  #/rad 

'W1 

Dimensional 
due  to 

-force  propeller  coefficient  #/rad^ 

X,.. 

Dimensional 
due  to  <£* 

/?£  -force  propeller  coefficient  #/rad^ 

or 

Kr  = 

Dimensional 
due  to  cf* 

^  -force  propeller  coefficient  #/rad 

= 

2 

Dimensional  if  -force  propeller  coefficient  #/rad 

due  to  A  cT*  / 

it 

Dimension.il 
due  to 

y*  -force  propeller  coefficient  #/ra.d 

-  Dimensional  J'  -force  propeller  coefficient  tf/rad 
due  to  A  r/'  rfj 
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Km 

Unit 

Dimensional  fa  -moment  propellor 
coefficient  duo  to  &cf 

0/  rad 

& 

ti 

Dimensional  -moment  propeller 

coefficient  due  to  c(4or  c/’j*’ 

0/ rad2 

'' 

Dimensional  -moment  propeller 

coefficient  due  to  (AcfJ' 

0/ rad2 

$ 

(I 

Dimensional  p  -moment,  propeller 
coefficient  duetto  eT, 

0/  rad 

*4  ■ 

Dimensional  <f  -moment  propeller 
coefficient  due  to  eQ, 

0/  rad 

^•4  ’ 

Dimensional  ^  -moment  propeller 
coefficient  due*  to  Z </•  oju 

0/  rad2 

A U  - 

Dimensional  j-  -moment  propeller 
coefficient  due  to  «£, 

0/rad 

N<,  ’ 

Dimensional  J  -moment  propeller 
coefficient  due  to  cf, 

0/  rad 

NaS-cf,  * 

Dimensional  p  -moment  propeller 
coefficient  duetto  &<f -  <ft 

0/rad2 

Propeller  thrust  parameter 

nondimensional 

Propeller  torque  parameter 

nondimensional 

V 

Propeller  power  parameter 

nondimensional 

P  = 

Power 

ft-#  /  sec 

^3 

11 

Power 

hp 

P  , 

Nondimensional  power 

r , 

Propeller  thrust  coefficient 

nondimen  sional 

/’  - 

propeller  torque  coefficient 

nondimensional 

Jr  ’ 

tan  IT 

nondimensional 

Zero  .subscript  denotes  trim  conditions, 
straight  and  level  flight 

j~  ami  sulmc  ript.s  ileuote  forward  and 

.all  conditions 

C,  ■•'ub.'ii  ript  denotes  propeller  rout  nil  In’ins 
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IV 

POSTULATED  TANDEM  PROPELLER  SUBMARINE  CONFIGURATION 


4.1  HULL 

Tin?  submarine  configuration  postulated  for  study  in  the  Tandem  Propeller 
Submarine  program  is  shown  in  profile  in  Figure 4-L  The  outline  was  taken 
from  Reference  fl)  and  is  purported  to  be  that  of  the  Thresher  submarine, 
SS(N)593.  However,  except  for  the  general  outline,  all  submarine  con¬ 
figuration  parameters  are  hypothetical 

In  ordci  to  simplify  the  stability  and  control  analyses  two  important  assump¬ 
tions  were  made  relating  to  the  submarine  hydrostatics  and  its  mass/inertia 
properties.  These  are: 

(1)  The  submarine  is  neutrally  buoyant  and  trimmed  at  ?ero 
angles  of  attack  and  sideslip,  in  straight  and  level  flight. 

The  centc  r-of-buoyancy  is  on  the  ^-axis  ^  feet  above 
the  cente r-of-gravity  (the  origin  of  coordinates). 

(2)  Mass  and  inertia  properties  are  calculated  on  the  basis 
of  a  homogeneous  prolate  spheroid.  The  sail  is  ignored 

in  mass/inertia  computations  but.  included,  where  possible, 
in  the  determination  of  hydrodynamic,  properties. 
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Thu  motacentric  height  ^  wa»  computed  on  the  basil  of  a  linear  interpolation 
between  submarines  of  tower  and  higher  displacement  (the  Albacore  and  the 
George  Washington).  The  distance  from  the  forepoint  to  the  c.g.  was  deter¬ 
mined  in  the  same  way.  The  propeller  planes  wore  located  as  far  from  the. 
c.g.  as  practical,  at  hull  stations  fore  and  aft,  where  the  diameter  is  approxi¬ 
mately  20  feet.  For  simplicity  the  distances an  in  Figure  1  were 

made  equal. 

The  physical  properties  of  the  postulated  TPS  hull  are  summarized  below: 

Symbol 


Submerged  displacement  :  4300  long  tons  (-  weight) 

Length  ;  275  feet  jt 

Max.  diameter  :  32  feet  d 

Fineness  ratio  :  8.6  £ 

Distance,  forepoint  to  c.g.  :  125  feet  ^ 

Distance,  prop,  plane  to  c.g.:  110  feet,  fore  and  aft  ^z 

Metacentric  height  :  1.0  feet 


4.2  PROPELLERS 

A  hub  diameter  of  20  feet  was  assumed  for  both  the  forward  and  aft  propellers 
and  each  propeller  was  assumed  to  be  equipped  with  16  blades  with  span  2  feet 
and  chord  1 .  5  feet. 

Propeller  rotation  is  clockwise  looking  forward  for  the  forward  propeller  and 
counterclockwise  looking  forward  for  the  aft  propeller.  Propeller  geometry 
is  defined  in  Fig.  4-2  and  the  important  propeller  physical  parameters  are 
listed  below. 
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20  feet 
24  feet 
16 


Symbol 


Hub  diameter 
Tip  diameter 
Number  of  Bladen 
Max.  rpm 
Blade  area 
Average  radius 
Average  diameter 


50 
3  ft2 
11  feet 
22  feet 


A/ 


Collective  and  cyclic  components  of  blade  pitch  are  shown  in  Figure  4-2. 
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V 

HYDRODYNAMIC  ANALYSIS 


5.1  GENERAL  REMARKS 

5.1.1  Task  Objectives 

The  hydrodynamic  analysis  (summarized  in  this  section)  has  the  three-fold 
objective  of  determining  the  (l)  stability  coefficients,  (Z)  control  coefficients 
and  (3)  added-mass  coefficients  of  the  tandem-propeller  submarine.  Each 
of  the  three  sets  of  coefficients  derive  from  the  force  interaction  between  the 
proposed  vehicle  and  its  surrounding  fluid  medium.  Traditionally,  however, 
the  stability  coefficients  are  defined  as  those  terms  that  yield  the  six-degree- 
of-freedom  forces  and  moments  acting  on  the  vehicle  as  a  result  of  the  total 
velocity  of  the  vehicle  relative  to  the  fluid.  The  control  coefficients,  on  the 
other  hand,  define  the  six-degree -of-freedom  forces  and  moments  produced 
by  whatever  control  mechanism  is  employed.  In  the  case  of  the  proposed 
vehicle,  the  basic  control  variables  or  inputs  are  the  collective  and  cyclic 
pitch  of  the  forward  and  aft  propellers,  respectively.  The  added-mass 
coefficients  yield  the  forces  and  moments  (for  all  six  degrees  of  freedom) 
that  are  caused  when  the  vehicle  is  accelerating  with  respect  to  its  surround¬ 
ing  Quid  medium.  Since  these  forces  are  proportional  to  vehicle  acceleration, 
both  linear  and  angular,  they  sum  with  the  inertia  reaction  terms  in  the 
equations  of  motion;  hence,  the  origin  of  tin?  term  "added  mass". 
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5.1,2  Problem  Areas 


The  most  straightforward  procedure  for  obtaining  the  hydrodynamic  charac¬ 
teristics  of  the  tandem-propeller  submersible,  (the  anumed  geometry  being 
given  in  Figure  4-1)  would  be  model  tests.  Lacking  test  data  such  as  would 
be  obtained  from  a  water  tunnel,  towing  tank,  or  whirling  arm  facility,  it  is 
necessary  to  resort  to  theoretical  predictions  baaed  on  fluid  mechanics 
analyses  and  to  extrapolation  of  empirical  data  obtained  for  similar  geometric 
forms . 

Since  the  tandem-propeller  submarine  concept  envisages  the  elimination  of 
the  usual  fixed  and  movable  planes  and  fins  at  the  aft  end  of  the  vehicle,  the 
hydrodynamic  analysts  resolves  to  a  task  of  predicting  the  force  and  moment 
properties  of  a  streamlined  body  of  revolution  as  modified  by  the  presence  of 
a  vertical  fairwater  and  two  counter-rotating  propellers  wrapped  circum¬ 
ferentially  around  the  forward  and  aft  ends  of  the  hull  (see  Figure  4-1).  A 
sizeable  number  of  practical  and  theoretical  problems  immediately  arise. 
First  of  all,  the  stability  coefficients  of  a  body  of  revolution  are  determined 
in  practice  by  those  properties  of  a  real  fluid  that  make  the  application  of 
potential -flow  theory  invalid.  The  flow  about  a  slender  body  of  revolution  is 
characterized  by  a  large  boundary  layer  and  flow  separation,  plus  a  wake. 
Thus,  it  is  not  possible. to  theoretically  predict  stability  coefficients  any  more 
reliably  than  one  can  predict  the  drag  coefficient.  Lacking  test  data  for  the 
specific  geometry  assumed  here,  it  is,  however,  possible  to  obtain  stability 
coefficients  by  extrapolating  from  the  body  of  test  results  obtained  by 
previous  investigators. 

In  addition  to  the  problem  of  interaction  between  the  flow  produced  by  the 
fairwater  and  the  three  dimensional  hull,  there  is  the  more  serious  problem 
of  hull-propeller  interaction.  From  Figure  4-1,  it  is  seen  that  the  flow  over 
the  hull  will  be  influenced  by  the  propellers  and,  concurrently,  the  flow  field 
at  the  propellers  will  be  influenced  by  the  presence  of  the  hull.  The  form 
and  magnitude  of  these  interactions  depend,  to  a  very  large  extent,  on  the 
magnitude  of  the  forward  velocity  ( (4  )  of  the  submarine. 
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For  tho  c.tt*«  of  high  forward  speed,  the  ;» h sumption  Is  made  that  the 

velocities  Induced  by  the  propeller  (both  forward  and  aft  of  each  propeller 

plane)  are  small  relative  to  the  forward  velocity,  .  Note  that  the  axial 

|  component  of  these  propeller-induced  velocities  will  Increase  the  axial  flow 

j  velocity  over  the  hull  above  that  due  to  tho  free  stream  velocity.  It  was 

estimated  that  tho  propeller-induced  axial  velocity  could  become  as  largo  as 

t  twenty  (20)  percent  of  tho  free-stream  velocity.  Since  the  existing  Reynold's 

.  7 

number  (approximately  10  )  indicates  fully  developed  turbulent  flow  (such  that 

the  rate  of  change  of  force  coefficient  with  velocity  is  small)  it  is  believed 

!  that  the  error  introduced  by  neglecting  the  influence  of  propeller-induced  axial 

and  tangential  velocities  on  hull  hydrodynamic  forces  is  small.  A  similar 

l 

i  conclusion  can  not  be  made  for  the  case  of  low  forward  speed  and  the  hover 

!  condition. 

In  addition  to  the  influence  of  the  propeller  on  the  flow  about  the  hull,  there  is 
the  influence  of  the  hull  on  the  flow  seen  by  the  propeller.  For  motion  of  the 
!  submarine  along  the  -JL  axis  only,  the  flow  field  is  symmetrical  at  the  plane  of 

i  the  propeller,  with  the  flow  velocities  in  the  inviscid  flow  field  near  the  hull 

being  higher  than  the  free  stream  velocity,  LC  •  This  increase  in  the  inviscid 
j  flow  velocity  is  also  accompanied  by  the  presence  of  a  boundary  layer  that  reduces 

the  fluid  velocities  seen  by  the  propeller  blades  near  the  boundary  of  the  hull.  In 
the  main,  the  departure  from  free-stream  conditions  is  small  for  axial  flow  per¬ 
mitting  the  justifiable  assumption  that  axial  motion  of  the  hull  has  no  significant 
influence  on  the  flow  field  existing  at  the  plane  of  the  propellers.  A  similar 
assumption  can  not  be  made  for  the  case  of  transverse  motion  of  the  hull  and 
the  influence  of  the  flow  field  (produced  by  transverse  velocities)  on  the  stability 
coefficients  of  the  propellers  is  discussed  below  in  Section  5.  2.  1.  2.  2. 

The  phenomenon  of  propeller  interaction  is  also  deserving  of  attention.  In  the 
case  of  high  forward  speed,  tho  previous  assumption  that  the  propeller-induced 
velocities  are  small  relative  to  the  free-stream  velocity  also  says  that  the  wake 
of  the  forward  propeller  has  little  influence  on  the  flow  field  existing  at  the  plane 
of  the  aft  propeller.  As  the  tandem  propeller  submarine  proceeds  trom  high-  to 
low-forward  speed  operation,  this  assumption  of  <ero  interaction  between  pro- 
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pullers  bucotnni)  untenable.  A  qualitative  discussion  of  the  influence  of  thu 
wake  of  Ihe  forward  propeller  upon  the  aft  propeller  is  given  in  Reference  (2). 
Some  of  the  propeller  interaction  phenomena  existing  at  zero  (hover)  speed  are 
discussed  below  in  Section  5.  2.  Z.  Z. 

5.2  STABILITY  COEFFICIENTS  OF  THE  TPS 

The  stability  coefficients  that  define  the  forces  and  moments  acting  on  the  TPS 
as  a  result  of  its  total  linear  and  angular  velocity  are  derived  separately  for 
the  cases  of  high  and  low  (or  zero)  forward  speeds.  It  will  be  seen  that  for  the 
low  speed  or  hover  case,  the  stability  characteristics  of  the  propellers  can  not 
be  derived  in  coctficient  form  nor  can  the  stability  forcos  and  moments  be 
separated  from  the  forces  and  moments  produced  by  the  control  variables, 
collective  and  cyclic  pitch. 

5.2.1  High-Speed  Case 

5.  2.1. 1  Hull  Plus  Fairwater  Contribution 

5 .  2 . 1 . 1 . 1  Evaluation  Technique 

The  drag  coefficient  of  the  TPS  hull  was  derived  from  the  experimental  data 
presented  in  Reference  (3).  This  reference  presents  data  for  a  series  of  bodies 
of  revolution  that  are  similar  to  submarine  hulls.  In  view  of  the  uncertainty  in 
exact  geometry  of  the  TPS  hull  form,  the  drag  coefficient  was  assumed  to  be  the 
average  of  the  drag  coefficients  presented  in  Reference  (3)  for  zero  angle  of 
attack  .  The  influence  of  speed  on  drag  was  not  determined  in  these  tests 
which  wore  conducted  at  a  constant  Reynolds  number  of  3.1  x  10^\.  Since  this 
Reynolds  number  is  sufficiently  large  to  insure  fully  developed  turbulent  flow 
(similar  to  that  over  the  full-scale  hull),  the  variation  of  drag  coefficient  with 


It  is  recognized  lliat  the  TPS  porf  jrmance  computed  for  a  specified  installed 
horsepower  will  he  a  sensitive  function  of  this  assumed  drag  coefficient. 
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velocity  can  bo  justifiably  neglected  in  tho  high-speed  case. 

As  mentioned  above  in  Section  5.  L2,it  is  not  practical  to  theoretically  evaluate 
tho  stability  coefficients  of  a  slender  body  of  revolution.  Hence,  experimental 
data  were  utilized  to  the  extent  that  such  data  were  available.  Reference  (4) 
presents  some  "build-up"  data  for  an  early  Fleet  Ballistic  Missile  (FBM) 
Submarine  configuration  with  a  fineness  ratio  (length  of  hull/maximum 
diameter)  of  11.  3.  The  hull  form  and  finonoss  ratio  were  deemed  sufficiently 
close  to  the  postulated  TPS  geometry  (fineness  ratio  of  8.6)  to  justify  the  use 
of  these  data.  Tho  stability  coefficients  of  the  hull  were  obtained  by  operating 
on  the  data  using  graphical  methods.  In  addition,  the  force  and  moment  con¬ 
tributions  of  the  assumed  fairwater  (see  Figure  4-1)  were  determined  using 
tho  low-aspect-ratio  hydrofoil  theory  presented  in  Refcronce  (5).  The  fair- 
water  (sail)  was  assumed  to  be  a  flat  plate  and  the  aspect  ratio  was  determined 
by  correcting  the  span  for  the  effects  of  hull  proximity.  The  computational 
details  are  not  included  in  this  report. 

5.  2. 1.1.2  Tabulation  of  Results 

The  results  of  the  above  analysis  were  placed  in  non-dimensional  coefficient 
form  in  order  to  facilitate  comparison  of  the  hydrodynamic  characteristics  of 
the  TPS  with  those  of  conventional  submarines  possessing  stabilizing  planes 
at  the  aft  end.  Table  5-1  summarizes  the  non-dimensional  stability  coefficients 
r  the  TPS  minus  the  fore  and  aft  propellers  and  compares  them  with  the  co¬ 
efficients  of  the  Albacore.  The  most  notable  difference  between  the  Albacore 
and  the  TPS  is  the  Albacore's  significantly  larger  damping  in  pitch  and  yaw 
(  //),£  and  \r'  )  and  vertical  force  (C  )  due  to  vertical  velocity  {^-  )  and 
pitching  velocity  (  ). 
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TABLE  5-1 


Non-Dimensional  Hull  Derivatives  (TPS  feAlbacore 
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TABLE  5-1  -  Tabulation  of  Nondiincnsional  Hull  Derivatives 
(TPS  ...  .d  Albacoro) 
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5,2. 1.2  Propeller  Contribution 

Any  departure  of  the  submarine  from  a  steady  course  (represented  by  a  fixed 
value  of  axial  velocity,  )  results  in  changes  of  flow  through  the  tandom  pro* 
pellers  and  consequent  force  and  moment  inputs  from  the  propeller  to  the 
hull.  Specifically,  when  the  propellor  has  motion  other  than  along  its  sp  axis, 
it  is  seen  that  the  advance  ratio  will  vary  around  the  circumference  of  the 
propeller.  These  variations  in  advance  ratio  will  produce  force  distributions 
on  the  propeller  disc  similar  to  the  distributions  produced  by  cyclic  pitch. 

Forces  and  moments  will  result  that  oppose  the  disturbance  motions  and 
various  cross-coupling  forces  and  moments  will  appear,  e.g.,  pitching  moment 
due  to  sideslip.  Although  an  analysis  has  been  performed  to  derive  the  stability 
coefficients  of  a  propeller  in  such  a  manner  that  the  analysis  is  valid  for  the 
complete  speed  regime  (see  Section  5.  2.  2),  the  work  of  Glauert  has  been  employed 
to  obtain  the  high  speed  stability  coefficients  of  the  large-diameter-to-hub-ratio 
propellers  employed  on  the  TPS. 

5.2. 1.2.  I  Results  Obtained  by  Glauert 

In  Reference  (6),  Glauert  derived  the  stability  coefficients  of  an  isolated  pro¬ 
peller  (no  hub  or  center  body)  for  the  case  where  the  disturbance  velocities 
(H  ,  ip  ,  jZr  ,  Cp  ,  &  ,  P  )  are  small  relative  to  the  free-stream  equilibrium 
velocity,  .  This  analysis  assumed  the  propeller  blade  angle,  £  ,  to  be  fixed 
and  the  load  distribution  on  the  blade  to  be  known.  Two  forms  of  engine -power 
function  were  considered.  One  function  states  that  the  ratio  of  engine  power  to 
engine  rpm  is  fixed  (i.e.,  P  / ZV"*!  ~  a  constant);  the  second  function  states 

that  engine  power  is  varied  to  maintain  propeller  rpm  fixed.  In  the  latter  case, 
the  speed-torque  characteristics  of  the  engine  do  not  enter  into  the  derivation  of 
the  stability  coefficients. 

Table  5-2  summari/.es  the  stability  derivatives  of  a  propeller  as  obtained  by 
Glauert  who  used  blade-element  theory  and  neglected  the  influence  of  inflow  and 
swirl.  Tile  lambda  parameters  shown  in  this  table  are  defined  by  Glauert  as  - 
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/V‘£  A  $f 

where  $  =  -  £4^  y 

and  T"  and/*  are  the  thrust  and  torque  coefficients,  respectively,  of  the  pro¬ 
peller.  Or  using  blade-element  theory,  these  thrust  and  torque  coefficients  can 
be  expressed  as  _ 

Since  propeller  rpm  will  be  considered  as  fixed  in  the  subsequent  stability  and 
control  analyses,  it  is  only  necessary  to  evaluate  the  A^and  X ^  parameters  by 
properly  formulating  the  derivatives  of  the  thrust-  and  torque -coefficient  ex¬ 
pressions.  Accordingly,  the  /\^  and  parameters  have  been  evaluated  as  - 

=  '  l  G.rSfrr)-n-GL[/-2S']+€Gto-Ypi-(s-Y)\ 

^  $  (s- rt-lc,,-  h o*(s-*y 

A0  =  /  f  3*  ,  j  1 

9  z  [F?  iti*  i$,(s-Y)  +  clll)+flGM*(s-*)*  j 

Note  also  that  the  parameter,  .  appearing  in  Table  5-1,  i»  obtained  from  the 
following  relationship  appearing  in  Reference  6: 
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1  Sj, 


blade-tip  radlua 

A/ for*  Jr  =  4  4fis£ 

blade- root  radlua  r 


blade-tip  radlua 

where  £  -NJ^^c^T 

blade-root  radiua 


<3  -  thrust  per  u  lit  length  of  blade 
T  -  radial  distance  along  blade 

In  order  to  apply  Glauert's  results  to  the  calculation  of  propeller  stability  co¬ 
efficients,  one  must  establish,  first,  the  trim  conditions  of  speed  (  tf0  ),  thrust 
and  torque  (  and  /^0),  blade  pitch  anulo  (  ),  and  blade  flight-path  angle 

(  'fp  )•  Also  required  are  the  slope  of  the  lift  curve  and  drag  coefficient  of  the 
individual  blade  elements.  These  basic  hydrodynamic  quantities  were  evalu¬ 
ated  by  assuming  each  propeller  blade  to  be  isolated  and,  effectively,  a  semi¬ 
wing,  The  curvature  of  the  wall  (hull)  was  neglected  and  the  lift-curve  slope 
was  evaluated  for  the  blade  geometry  described  earlier  by  making  the  classical 
aspect  ratio  correction  (see  Reference  (7)).  The  following  blade  hydrodynamic 
characteristics  were  obtained: 

Q  --  3. 59 

Cn  s  .015 

uo 

f,  =  .119 

5.2.  1.2.  2  Effect  of  Transverse  Flow  Around  Hull 

As  mentioned  in  Section  5.1.2,  the  flow  field  at  the  propeller  plane  when  the  pro¬ 
peller  has  a  transverse  velocity  relative  to  the  fluid  will  be  altered  by  the  pre¬ 
sence  of  a  body  of  revolution  centered  within  the  propeller  disc.  The  local  flow, 
as  caused  by  the  presence  of  the  center  body,  can  be  approximated  reasonably 
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well  by  evaluating  the  two-dimensional  potential  flow  around  the  circular  cross 
section  of  the  hull.  The  local  flow  velocity  caused  by  either  a  heaving  or  side* 
slipping  velocity,  /at-  or  AT  ,  is  given  by  (see  Reference  (8))  - 

-"''local  ‘Z*si,s'»tr 
^  local  •'  2  «" 

where  and  rsfi  refer  to  the  velocity  of  tho  propeller  as  defined  at  the  origin 
of  an  axis  system  fixed  in  the  plane  of  the  propeller.  Note  that  the  local 
velocity  varies  around  the  hull  as  a  sine  or  cosine  function  of  the  blade  azimuth 
angle,  O”  ,  in  the  same  manner  that  the  free -stream  velocities,  and  /Vp  , 
alter  the  tangential  velocity  in  a  blade  element  when  no  center  body  is  present. 
In  view  of  the  identical  functional  form,  the  integration  process  involved  in 
summing  the  forces  on  one  blade  over  a  complete  traversal  of  azimuth  angle, 

<T  ,  is  not  altered  but  there  is  a  twofold  increase  in  stability  forces  caused  by 
transverse  velocity  when  a  center  body  is  present.  This  correction  has  been 
applied  to  the  derivatives  •  and  shown  in  Table  5-2. 


5.2.  1.2.3  Tabulation  of  Results 


The  above  presented  stability  coefficients  apply  to  one  propeller  only  and  to 
propeller  axes,  where  these  axes  are  parallel  to  body  axes  but  have  their  origin 
at  -96  =  ,/->  ,  fore  and  aft.  Thus,  the  following  relationships  must  be  used  to 
sum  the  contributions  from  the  fore  and  aft  propellers  such  that  the  combined 
coefficients  or  derivatives  refer  to  motions  of  the  c.g.  of  the  TPS: 

<Tp]  =  ML  ;  ^]=  &  +  ^4  ;  rZ'  1  =  nr- 

-  "f 
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Hr'1 :  e4/6 


X  J 
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*1  =  7 
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TABLE  5-2 

Tabulation  of  Propeller  Derivative*  (One  Propeller) 

Sign  Change  for 
Counterclockwise 

Clockwise  Rotation  Looking  Forward 

Rotation 

Xu  m  +2  ^Vofi  (}{  )  /u.c 

No 

x?  *+2Fx.f i(v/n. 

Yes 

Yv  -  -2lM*oki 

No 

Vj  *+KJj,/u0 

z£  =  r^ 

Yes 

/ 

K  -  t? 

/ 

No 

Yes 

^lv  *  “2.  *Rt 

Yes 

3  fei.  R<  t>  D*  *^r/(lo 

No 

H«£  -  ai^- 

1 

/ 

f  ('TO 

CONSTANT 

ComstAMT-A. 

f,00 

r^l-f  c\p)- / 1 4- "Rq 

MnV 

T^l-  Hr)  /•+ 

\  —  Rt 

fjfh) 

7\p6  -  /l-b  ftp-  R «, 

1“  rl* 

VcO 

/ l-hTVp-*  ^ 

TABLK  5-2  -  Tabulation  of  Propeller  Derivatives  (One  Propeller) 
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The  total  inability  forces  and  moments  produced  by  the  propellors  are  thus: 


I  %  =  £[(*++4)+  xj](*-  K)  +  y/l/ 

IzP=  7  v 

Jkp  «  Kf]  rP  +  K r ]  * 

p  f-y.*  N  -if,i 

I/i,  =  ^  *  KJ("  -  v* ) 4 

+ + ®4  >4  vA  +  <]/4  +• 

Z^  =  /O + ^  -  < >  -  M  )4  +  < ](■=-  -  &A) 

+  )+y/J/4  -  y/>4  +  <£/ 


U  it  is  auumed  that  the  fore  and  aft  propeller*  are  symmetrical  and  are 
operated  at  a  constant  propeller  speed  with  the  forward  propeller  turning  clock¬ 
wise  and  the  aft  propeller  turning  counterclockwise,  then  insertion  of  the  proper 
propeller  coefficients  (from  Tablo  5-2)  into  the  above  force  and  moment  sum¬ 
mations  produces  the  following  results: 
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r  4lvj(/-\a) 

L  JT7F~ 

¥ 

II 

K! 

r  4 / u* to- xj 

L  _n.Rz 

■]' 

M 

l^sl 

r  4 IM^I(I-Xa) 

L  ji 

-]* 

II 

l/\] 

.  r  a &£!.(/ 

L  ji 

+[2  k*  ^  ~ 

^o-\)ey 

4-HJO-  \*)Ak 

r lg* 


10 


AO- It.  l-l  ■  V- 1 


The  above  force  and  moment  summations  are  expressed  in  terms  of  the  trim 
thrust  and  torque  (/ and/^)  of  a  single  propeller  and  in  terms  of  other  trim 
quantities  such  as  H0  •  J~L  ,  A  r  ,  »  and  •  For  unbalanced  operation 

of  the  propellers,  the  above  summations  are  not  valid  but  must  he  redorived 
taking  into  account  differences  (between  the  forward  and  aft  propellers)  in 

,  jO  ,  and  the  lambda  parameters  which,  in  turn,  are  functions  of 
y  and  (f  • 

5.2. 1.3  Combined  (Hull-Propeller)  Coefficients 

The  stability  characteristics  of  the  hull-propeller  combination  can  be  put  into 
derivative  form  such  that  the  equations  of  motion  reduce  to  their  simplest  form. 
Since  the  hull  stability  derivatives  have  been  tabulated  in  non-dimensional  form 
and  the  propeller  forces  and  moments  have  been  evaluated  in  dimensional  form, 
the  decision  was  made  to  dimensionallze  the  hull  derivatives  and  to  sum  these 
derivatives  with  their  propeller  counterparts.  Accordingly,  a  set  of  dimensional 
stability  coefficients  is  obtained  based  on  perturbation  velocities  such  that  the 
total  hydrodynamic  force  and  moment  acting  on  a  TPS  (exclusive  of  control  forces) 
can  be  expressed  as  follows: 

Xx  =  +  XK  H 

ly=  y,*  +  yrv  + 

IZ-Z^  +  Zpe 
Y_K  -  +  Kf?  +  Krf+  /%«**"** 

J_V\  =  /%/*  +■  f  Sl"  * 

XW*  Nj?  +  A lf$  +  /Vfe  +  Nr  Y 

The  derivatives  in  the  above  equations  are  defined  and  tabulated  in  Table  6-2. 
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5.2.2  Hovering  Casa 


5.2.2.  1  General  Remarks 

The  propeller  forces  and  moments  plus  hull  forces  and  moments  produced  as  a 
function  of  velocity  disturbance  variable:)  for  the  case  where  i(0  is  small  or 
zero  are  derived  in  this  section.  It  should  be  noted  that  the  below  hydrodynamic 
analysis  is  preliminary  and  should  be  considered  tentative  until  such  time  that 
these  results  can  be  checked  further.  One  drawback  with  the  analysis,  at  the 
present  time,  is  that  the  definitions  and  terminology  used  in  Reference  (9)  are 
employed  rather  than  that  employed  in  the  remainder  of  this  report.  Neverthe- 
less,  the  results  serve  to  indicate  the  general  nature  of  the  problems  that  will 
be  encountered  in  performing  a  stability  and  control  analysis  for  the  hover 
condition. 

The  analysis  made  to  evaluate  the  low-speed  propeller  forces  is  based  on  an 
analytical  procedure  believed  to  be  sufficiently  general  that  it  permits  the 
derivation  of  propeller  forces  and  moments  throughout  the  speed  regime.  Only 
the  low-speed  solution  is  given  here,  however.  To  avoid  mathematical  detail, 
only  th<>  bare  outline  of  the  procedure  is  given  with  the  final  result  obtained  for 
the  X  -force  presented  in  the  body  of  the  report.  The  remaining  five  force  and 
moment  expressions  are  given  in  an  appendix. 

A  tentative  formulation  is  made  of  the  forces  and  moments  produced  by  the  hull 
at  low  speeds.  When  this  analysis  is  refined  and  reduced  to  a  quantitative  basis, 
it  is  anticipated  that  these  body  (hull)  forces  will  be  small  relative  to  the  forces 
produced  by  the  tandem  propellers. 

5. 2. 2.  2  Propeller  Forces 

With  reference  to  the  following  sketch,  blade -element  theory  yields:. 
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D(r)--i  PA  Vi*)  [cDc  +  f  c  i  SM  -  tn'  y£]J‘ 


where  =  total  (local)  water  velocity  relative  to  the  blade, 

ft/  sec. 

\Z(t)  -  local  tangential  velocity  relative  to  the  blade  in 
the  plane  of  the  propeller,  ft/ sec. 

}(*)  -  local  axial  velocity  relative  to  the  blade,  ft/ sec. 

In  contrast  to  the  assumptions  made  for  the  high-speed  case,  namely,  that  (1) 
propeller  forces  caused  by  disturbance  velocities  are  independent  of  the  blade 
angles  of  attack  caused  by  changes  in  blade  pitch  angle,  S  ,  and  (2)  changes  in 
blade  angle  of  attack  produced  by  pitch-angle  change  cause  variations  in  pro¬ 
peller-induced  flow  that  are  negligible  relative  to  the  free-stream  velocity,  l(0  , 
it  is  necessary  to  include  all  sources  whereby  ,  V  •  ^  •  ar>d  ®  vary 

around  the  circumference  of  the  propeller  as  a  function  of  the  blade  azimuth 
angle,  O'*  .  Both  cyclic  pitch  and  body  (hull)  motions  will  cause  perturbations 

in  £  and  V  and,  hence,  in  the  total  velocity,  V  .  Cor.jrtferation  of  the 

prevailing  kinematics  and  geometry  yields  the  following  for  the  local  axial  and 

>x 

tangential  velocity  components: 


>;« 

The  zero  reference  for  the  blade  azimuth  angle,  <J~  ,  is  the  positive  dt^axis 
as  in  Reference  9.  ^ 
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<4  +  siftr  -rAcos<r  {5. 

V^)*  (jTL-^R  +  (/V*  +  r'fz  )St'*  *  -  (.*/-- )  r«s  °“ 


where 


=  propeller-induced  axial  inflow  velocity  resulting 
from  momentum  chango  through  the  propeller 


Note  that 


V(<r)  =  ti?f 


(5-4) 


When  C(  =  t(0  (large  value)  and  /%r  -^W  =  -f>  =  = 1  X  =  0  '  .fc:  is 

assumed  negligible,  and  Equations  (5-3)  and  ( 5-4)  reduce  to 

#0  =  tiritt. 

Equations  (5-1)  and  (5-2)  for  lift  and  drag,  together  with  Equation  (5-3),  yield 

the  instantaneous  force  on  a  single  blade  as  a  function  of  blade  azimuth  angle, 

where  the  blade  angle  of  attack  and  the  total  water  velocity  relative  to  the  blade 

vary  continuously  around  the  circumference  of  the  hull.  Since  the  frequencies 

at  which  the  hull  responds  as  a  rigid  body  are  extremely  low  compared  to  the 

frequency  of  propeller  rotation,  it  is  evident  that  only  the  average  force  exerted 

>> 

by  a  propeller  blade  will  affect  the  hull  motions  .  The  average  force  per 
revolution  of  a  single  blade  times  the  number  of  blades  is  the  game  as  the 
summation  of  the  forces  on  each  blade  at  a  given  instant  of  time.  For  example, 
the  total  average  axial  force,  ,  can  be  obtained  by  the  following  integration: 

M  Cos  j  /(a-)]  -  DM  sin  f  efer  (5-5) 

o 


where 

N- 


number  of  blades 


Statement  is  true  in  the  stability  sense.  Hull  vibrations  are,  however,  very 
much  a  function  of  the  periodic  propeller  forces  even  though  the  mean  value 
may  reduce  to  r.ero. 
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The  integration  procedure  indicated  in  Equation  (5-5)  would  be  u»ed,  in  similar 
fashion,  to  derive  the  remaining  force  and  moment  component*  produced  by  the 
propeller. 

Examination  of  the  blade-force  equations  (lift  and  drag)  indicates  that  they  are 
nonlinear,  for  example,  in  ^  and  V  •  The  above  indicated  integration  is 
facilitated  by  linearizing  the  terms  in  the  force  equations  about  the  trim  speed. 
Expansion  of  (i.e.,  the  inverse  tangent  expression,  tan'  )  in  a  Taylor's 

series  and  linearizing  by  retaining  only  the  constant  and  linear  terms  leads  to: 

Un'  4  *  t*i' '(nR  ^  f V  -  jrR] 

o 

where  is  the  trim  equilibrium  speed.  At  the  trim  advance  ratio,  jfc  ft  > 
of  1.0  this  expansion  gives  good  accuracy  in  the  range  0.  85<  J*Q^1. 15.  Simi¬ 
larly,  the  total  velocity,  ,  can  be  approximated  by  means  of  the  binomial 
expansion: 

-  yJJTJif 


In  the  case  of  the  inverse  tangent  expansion,  the  quantity >3  treated 
as  a  small  quantity  and  in  the  case  of  the  total  velocity  expansion,  the  quantity 
iar-  fe§rfj  is  treated  as  a  small  quantity.  In  these  expansions,  only  first 
powers  in  the  motion  variables,  ,  /XT’ ,  an <\yu/“<  are  retained. 

Products  of  these  terms  are  also  neglected  in  the  derivation  of  the  final  force  and 
moment  expressions. 
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Evaluation  of  the  induced  inflow  velocity,  £/  ,  can  be  carried  out  by  meana  of  a 
combined 'lifting  line  and  propeller  momentum  analyeie.  Aa  a  conaequence, 
can  be  expreaaed  aa  a  function  of  the  trim  collective  pitch,  ,  where  the  total 
pitch  angle  of  the  blade  ia  defined  (similar  to  Reference  (9))  aa: 

S(<r  )  -  SQ  +3,  Si  ??  cr  f  cos  cr 

The  above  procedure  ia  equivalent  to  atating  that  the  induced  velocity  and  thruat 
change  with  the  collective  pitch  angle  without  lag.  Effecta  of  the  cyclic  pitch  on 
the  induced  velocity  are  neglected  and  the  caacade  effecta  are  neglectod.  Note 
that  the  adequacy  of  the  firat  aaaumption  (neglecting  the  effecta  of  cyclic  pitch) 
dependa  on  the  operating  condition.  It  is  a  good  aaaumption  for  high  apeeds  where 
the  propellers  are  operating  at  high  thrust  levels  and  is  a  reasonable  assumption 
for  that  low  speed  (or  hover)  case  in  which  the  propellers  are  operated  at  high 
thrust  levels  by  "bucking"  one  propeller  against  the  other.  It  is  possible  that 
further  analytical  work  could  permit  one  to  relax  these  assumptions. 

For  the  case  of  low  or  zero  forward  speed,  the  terms  in  the  blade  force  equations 
are  linearized  about  tj0  =  Q  .  The  expansions  and  integrations  have  been  carried 
to  completion  for  the  low  speed  case  and  the  results  for  the  six  forces  and  moments 
are  presented  in  an  appendix.  The  result  obtained  for  the  ^  force  is  given  below. 
Note  that  this  expression  is  not  valid  when  (  C{.  +  )  approaches  zero  -  a  possible 

operating  condition  when  the  thrusts  of  the  propellers  are  opposed  at  low  speeds. 


U. 
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whore 


_L=  L<i«  f-  /  + 

Rfi  8-2$,  i  1  + 


/Vc 

/? 


£  s  mean  blade  chord 


yielding 


Xl 

KJX 


0.858 $0  -  0.600  S* 


Examination  of  Equation  (5-6)  shows  that  there  is  coupling  between  control 
deflections  and  motions  of  the  hull.  As  will  also  be  shown  later  for  the  high 
speed  case,  there  are  cross  couplings  between  the  various  components  of 
control  deflections.  Accordingly,  Equation  ( 5-6),  (plus  the  Appendix)  indicates 
that  a  stability  and  control  analysis  of  the  hovering  case  will  be  extremely  com¬ 
plicated,  unless  the  indicated  coupling  effects  can  be  shown  to  be  small  for 
reasonable  values  of  control  inputs.  These  questions  will  be  resolved  during 
the  remainder  of  this  study. 


5.  2.  2.2.1  Possible  Modes  of  Operation  at  Low  Speed 

Considerable  flexibility  exists  in  the  choice  of  a  low  speed  mode  of  operation. 
These  include  at  least  the  following: 


(D 

<*) 

(3) 

(*») 

(5) 


Front  prop  thrusting,  rear  prop  non-rotating  but  trimmed 
to  counteract  torque 

Rear  prop  thrus.jjng,  front  prop  non  -  rotating  but  trimmed  to 
counteract  torque 
Both  props  thrusting 

Front  prop  thrusting,  rear  prop  counter-thrusting 
Rear  prop  thrusting,  front  prop  counter-thrusting 


These  cases  may  have  quite  different  stability  and  control  characteristics  and 
probably  are  quite  different  with  respect  to  vibrations  that  would  be  induced. 
This  latter  point  can  be  visualized  by  considering  two  possible  modes  of  vertical 
translation  at  zero  axial  speed  as  shown  on  the  following  sketch: 


In  Sketch  (4a)  the  wakes  of  the  propellers  approach  each  other.  This  would  cause 
a  confused  vortical  distribution  near  the  hull  with  possibilities  for  exciting  hull 
vibrations  and  propeller  vibrations. 

In  Sketch  (5a)  the  flow  over  the  major  portion  of  the  hull  is  non-vortical,  and  the 
wakes  from  the  two  propellers  do  not  approach  each  other.  Hence,  the  latter 
appears  to  be  the  more  desirable  from  the  vibration  viewpoint. 
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5. 2.  2.  3  Hull  forces  and  Moments 

The  hull  forces  and  moments  for  the  low  speed  or  hovering  cases  are  different 
from  those  for  the  hl^h  npeed  case.  In  the  low  speed  case  the  velocity 
\J  -\J  H*  +/tr*  f/cr  in  very  small.  The  lift  (1.  e. ,  side  force)  generated  by 
the  fairwater  will  be  extremely  nonlinear  and  is  omitted  here.  The  predominant 
terms  probably  arise  from  the  cross  flow  velocities.  /j*  and ></'*.  In  this  caso, 
viscous  cross  flow  theory  can  be  used  and  the  drag  components  resulting  from 
the  velocities  considered  independently.  These  forces  will  also  be  nonlinear 
and  small  but  at  least  they  can  be  defined  with  moderate  accuracy.  Their 
inclusion  is  necessary  for  the  definition  of  trim  conditions  when  n  -  O  .  The 
following  force  expressions  are  obtained: 

-  i  (-Dt PAtf'1  +  -k^t>3PA\ 

•  ~kcD,P 


where: 


P 


-  viscous  skin  friction  drag  coefficient  for  the  hull 
=  projected  frontal  area  of  hull  plus  sail 

=  velocity  component  along  *yL  axis 
=  density  of  fluid  medium 

=  viscous  cross  flow  drag  coefficient  for  the  hull 

=  projected  side  area  of  the  hull  (assumed  to  be  a  body 
of  revolution) 

=  drag  coefficient  for  the  sail  in  transverse  flow  (assumed 
to  be  a  flat  plate) 

-  projected  side  area  of  the  sail. 
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.iHiC.y/3 1* 


r  I  >J4'C,a*iS  S*1 


Then 


A,  *  ILjjz. 

Az  ''ft Mb 

+• o 

A,*J-c 

3  maximum  hull  diameter 
J&*  -  height  of  tail  above  hull 
t  -  thickness  of  sail  airfoil 
d(*)  3  diameter  of  hull  at  station  rf. 
C  =  chord  of  the  sail 

+s&t) 

^hviI  ^  Pss’JdMJ?  +  -fa 

z*%;\  ic*k  (+)<** 


and  the  terms  like  *  Xir  '  '  '  aro  zero* 


The  moment  expressions  become 


K  - 

n  - 

/V  ■ 

where  ^  = 


Rf  Z  ^ 

~k  C,  fa*  ft<K+)  <tV 

i  \  +R&  i  CDj  Pr/J-C 

normal  distance  from  axis  to  confer  of 


pressure  of  the  sail 


=  normal  distance  from 
pressure  of  the  sail. 


axis  to  cento r  of 


Af«  -  It.  1-1  •  V  - 1 


-s-jf  ,  ji^s-  =<• 


Thu  drag  coefficients  must  be  properly  evaluated,  taking  into  account  the  con¬ 
figuration  of  the  hull  and  the  Reynold's  number.  It  would  be  very  desirable  to 
have  experimental  data  for  the  drag  coefficients  but  in  view  of  the  lack  of  such 
values,  classical  drag  data  for  flat  plates  and  cylinders  must  be  used  .  Also 
the  expression  for  4  )  should  conform  to  the  desired  hull  form.  A  properly 

chosen  ellipsoid  of  revolution  would  probably  suffice. 

5.  3  CONTROL  COEFFICIENTS  OF  THE  TPS 

The  control  coefficients  of  the  tandem  propellers  are  derived  below  for  the 
case  where  the  TPS  operates  at  a  high  value  o *  forward  velocity,  tf0  .  (It 
should  be  recalled  that  the  control  forces  for  the  low  speed  or  hover  case  can 
not  be  expressed  independently  of  the  propeller  stability  forces.  Consequently, 
ho  analysis  of  the  low-speed  control  forces  is  given  in  this  section,  since  they 
were  treated  in  Section  S.  2.  2.  2  above,  albeit  in  a  preliminary  way. ) 

5.3.1  Evaluation  Technique 

The  propeller  forces  and  moments  arising  from  collective  and  cyclic  pitch  control 
have  been  derived  (using  blade-element  theory)  by  the  Netherlands  Ship  Model 
Basin  and  presented  in  Reference  (9).  The  derivation  is  briefly  reviewed  here 
in  order  to  obtain  expressions  that  are  consistent  with  the  terminology  adopted  by 
CAL  and  to  further  account  for  the  different  zero  reference  used  by  CAL  in  de¬ 
fining  the  propeller  blade  azimuth  angle,  G" 

Consider  the  hydrodynamic  forces  on  one  blade.  Using  blade  element  theory: 


In  the  absence  of  experimental  data,  one  should  base  Cpt  on  the  skin  friction 
drag  and  total  surface  area,  on  drag  data  for  cylinders,  and  Cp^  on  the 

drag  of  a  flat  plate  normal  to  stream  direction. 


11 


AG  - 1  (>  1-1  -  V  - 1 


D  •  irv*(<^+f!<£S)A 

where  ]/  ,  the  blade  velocity  with  respect  to  the  water  is 
given  by  (see  Figure  4-2): 

S/l= 

and  straight  and  level  (light  at  high-speed  is  assumed. 
Interaction  between  the  propellers  is  neglectod  as  are 
inflow  velocity  and  swirl. 


The  lift  and  drag  forces  are  first  resolved  in  directions  parallel  to  >  and  tangent 
to  the  average  blade  circle  (in  tho  propeller  plane)  and  then  along/*  •  y  •  'ty  to 
yield  expressions  for  the  propeller  control  forces.  Using  the  notation  and 
symbols  of  Figure  4-2  the  instantaneous  forces/moments,  written  in  propeller  axes, 
are  given  by: 


£ 

N 

4 

4 

r%. 


=  ~/+  I Frl  coS 


=  Vi  /frl 5<77cr 

=  "/t  I  Frl  ft 
My  =  “/-  /v  R  COS  cr 


(5-9) 

(5-10) 

(5-11) 

(5-12) 

1 

(5-1  1) 


My  ---/-  Ff  R  <r 

in  which  the  notation  ’/t  designates  tho  sign  of  the  term  according  to  whether 
the  propeller  in  question  is  forward/aft,  and  and  f  Rr  f  are  the  lift  and  drag 

resolved  along  ^4  and  tangent  to  the  average  radius  ft  .  As  in  tho  case  of  tin- 
stability  forces  in  5,2.  1.2,  3,  balanced  operation  is  assumed.  If  this  is  not  the 
case,  differences  between  the  forward  and  aft  propeller,  must  bo  taken  into 
account  in  (5-7)  thru  (5-13). 


ssaaa>a,swfcr!w*--y‘1'1  "y * 
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By  substituting  (5-7)  or  (5-8)  into  (5-9),  (5-10),  (5-U),  (5-12),  *0(1(3-13), 
letting  oC  Sin  O- -hfiiCosO* »  and  integrating  over  one  cycle  with 

respect  to  the  propeller  angle  0~,  the  equations  defining  the  average 
forces/ moments  can  be  written  in  body  axes  as; 

X:  Ac=  %  \  pVAN  j c^tbsvt,-  j j  \ 

Y*  Yc  =  "4  ^  f>  Am  i  sw v + rosif  i  o.‘  b.  j 

l-  ic  7+  |  f  Am  jC^ftnYf  (aiXf.fcJ  U  ?  Si 

(5-14) 

K:  Kc --  -/f  ^  VlAMR  [C,  W  t.  +G.sSffQ.  ffC*  j Cc  ii+ii  |  J  / 

M:  M ifVVRfq^x  -4  iju 

N:  Ne-y-ipAMR^l^y-SmY^C^U/S,  +/_  Yci^ 
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Some  interesting  properties  of  the  cyclicly  pitched  propeller  can  be  seen  from 
Equatlom(5-l4).  In  particular!  attention  is  called  to  the  following: 

(1) 


(2) 


(3)  The  forward  (clockwise)  propeller  produces  a  negative 
roll-moment  which,  for  balanced  operation,  is  cancelled 
by  a  positive  roll -moment  due  to  the  (counterclockwise) 
aft  propeller.  As  in  the  case  of  thrust,  there  is  a  coupled 
roll-moment  due  to  cyclic  pitch.  Notice  also,  that  since  the 
coupled  moment  is  due  to  induced  drag  it  is  independent  of 
the  signs  of  ot  and  5^  . 

(4)  The  pitch  and  yaw  moments  /%  and  contain  pure  couples, 
generated  at  the  propellers  (as  indicated  by  Equations  (5-12)  and 
(5-13)  and  due  to  the  fact  that  /y,  is  not  axial),  as  well  as  moments 
generated  by  the  propeller  forces  ZC  and  Yc  •  As  will  be  seen 
later  these  pure  couples  are  small  relative  to  the  terms 

and  • 

5.3,2  Reduction  to  Coefficient  Form 

The  control  forces  and  moments  given  by  the  set  of  Equations  (5-14)  can  be 
reduced  to  coefficient  form  to  facilitate  the  ensuing  stability  and  control  invest! 
gition.  The  eventual  numerical  evaluation  of  these  coefficients  will  demonstrat 
in  part,  the  extent  to  which  the  various  control,  cross  couplings  are  significant. 


Thrust  is  derived  primarily  from  collective  pitch,  but  is 
diminished  by  the  profile  drau  term  UATf  Cdq  and  the  induced 
drag  term  5 / 71  Vf,  The  latter  is  a  iunction  of 
collective  pitch  and  cyclic  pitch  8,  and  . 

y  force  is  obtained  by  applying  cosine  cyclic  pitch  (  S^) 
but  is  dependent  also  upon  the  level  of  collective  pitch 
(  69  ).  (A  similar  remark  holds  true  for  Z  force). 
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A  more  conclusive  answer  to  what  degree  these  cross-couplings  are  significant 
is  obtained  from  the  simulation  stud/  discussed  in  Section  VII. 

In  Equations  (5-14)  let  h0  3  80  +A$  -  Ya  ,  l.e.,  the  angle  of  attack  due  to 
collective  pitch  is  equal  to  the  trim  collective  pitch,  .%  ,  plus  the  change  in 
collective  pitch,  A  S  *  minus  the  trim  flight  path  angle  of  the  propeller  blado, 
^  .  Next,  perform  a  partial  differentiation  of  Equations  (5-14)  with  respect 
to  each  of  the  control  inputs  AS  ,  8,  ,  and  8z  •  The  method  is  illustrated 

below  by  differentiating  the  equation  for  the  X  control  force  and  defining  the 
resulting  coefficients.  (The  coefficients  resulting  from  differentiating  the  re¬ 
maining  five  equations  will  be  tabulated  without  discussion). 

On  letting  ~  $0  +  A$mTfo  •  basic  expression  for  thrust  of  the  forward 
propeller  is: 

Xc  =  i(’V’Atifo/‘S  Y'[S0+aS-X,l-si* 

Noting  that  becomes: 

Xc  =  cat  K[itAS-rJ  -  c0o-rA*  • 

[ ($0+Ai-xf+ 


Letting 


war — „ 

3  (A  r) 


=  Xa!  5  *  -OS) r  -  Xf 


UM  -  v 


...d  CX c-jL  -  X 

\  /rn  A 


At ;  - 1 «»  j  i  -  v  -i 


and  taking  the  partial  derivative*  indicated,  the^  -force  coefficients 
become  (after  the  substitution  o(Q  3  S0  -  ^  is  made): 


For  simplicity,  let 


The  force  and  moment  coefficients  for  the  remaining  axes  can  be  arrived  at  in 


the  same  way.  These  are  listed  below  (for  the  torward  propeller): 


V-  - Lev? AN 


■I  6 
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Using  the  above  force/moment  coefficients  the  following  set  of  control  force/ 
moment  equations  result.  Balanced  operation  is  assumed,  with  forward  pro¬ 
peller  and  aft  propeller  coefficients  being  equal  in  magnitude. 
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Yc  ~  tg  (S^~  Sif)  +  %f.f  (A&f 

Ze=Zg(S,f-*k)  +  Z&s-S  (&sf  sif- 4  sf  5'i  )  (5.l5) 

Kc=  Kas  (&%  -Mi  )  +  Ksi(%-%t%-%) 

/%  - 

.  /  /  .  ,  +t1ts.sA  (*$%  +**<*;  .) 

Nc  -  /vs (if Xg-s -4 

Notice  that  except  for  the  trim  thrust  {the  trim  moments  cancel  for  balanced 
operation)  the  above  set  of  equations  describe  porturbation  forces/moments  since, 
i i  &$  -  St  -  ~  0  •  they  reduce  to  zero. 

If  the  trim  term  in  5-lS)  is  ignored  for  the  moment,  the  control  forces  and  moments 
are  seen  to  consist  of  a  set  of  primary  control  terms  and  a  set  of  cross-coupled  or 
cross-axis  terms.  Let  the  control  inputs,  forward  and  aft,  be  defined  in  the 
following  way: 


/yL  -axis  (force)  control  input  = 

£ 

-  ASj. 

^-axis  (force)  control  input  = 

H 

-  $2 f  ~  &<? 

^  -axis  (force)  control  input  = 

r  -  c T,  3 

zjt  -axis  (moment)  control  input  = 

Sk 

=  A$  +  -A%Z 

^-axis  (moment)  control  input  = 

■z.  $lf  +  Sli 

^-axis  (moment)  control  input  = 

Sn 

*  <% 

The  primary  control  terms  are  then: 
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r 


X*  ^ 

y*  y$ 

K : 

/y.#  ms  *« 

H*  h'f  5y 


(5-17) 


and  the  undesirable  cross-coupling  terms  are: 

X:  i  V  W +%)  +  *  XrfS*  +  s*  +  #  *  # ) 

y*  ^  ) 

^  ^s-s*  (X  ^  *  ^/c  ) 

/<;  a  fH )  * Kf* (sff«  +w  (5’l8) 

xy»  ~ z t st*(fN) + ^Lg^cs^f^ + ^) 

/V«  ^  ^  ^  4<  f  $,  ^  ^ ^ ) 


Notice  that  most  of  the  coupling  terms  are  non-linear,  since  they  involve  the 
squares  and  products  of  the  control  inputs. 
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5.4  ADDED-MASS  COEFFICIENTS 

For  the  purposes  of  this  study,  the  moment*  of  inertia  of  the  TPS  were  computed 
by  assuming  the  aubmarine  to  be  a  homogeneous  ellipsoid  of  revolution.  The 
coefficient*  of  apparent  maaa  and  inertia  wore  derived  from  classical  sources  in 
the  literature,  again  assuming  the  submarine  to  be  a  prolate  ellipsoid.  The 
numerical  results  are  tabulated  in  Table  6-1,  whore  the  terminology  employed 
for  defining  the  real  and  "added"  masses  and  inertias  in  that  established  by 
SNAME. 
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VI 

EQUILIBRIUM  EQUATIONS  AND  NUMERICAL  DATA 


6.  1  THE  SIX-DECREE-OF-FREEDOM  EQUATIONS  OF  MOTION 
In  Section  V  the  propeller  and  hull  hydrodynamic  forces/momenta  were  written 
in  perturbation  form  using  dimensional  derivatives  which  incorporated  both 
the  hull  and  the  propeller  effects.  These  forces  and  moments  were  defined 
as:  M  •  The  propeller 

propulsion  and  control  forces/ moments  were  also  discussed  in  Section  V  and 
designated  Xc  .  Yc  .  2C  *  Kc  »  Me,  *  the  gyro¬ 

scopic  terms  are  added  to  these  two  groups  of  body  forces/moments  and  the 
results  are  equated  to  the  mass/ inertia  reaction  terms,  the  complete  six- 
deg ree-of-freedom  equations  of  motion,  in  abbreviated  form,  become: 


X*  1ri,U->vi4vu-+tV)J^W‘ •+•  X<t  ~'- 

V;  v^'Cr+in/u-TnjpvV-  £.f  +■  Yc 

Z:  TVjti-^j.U  +  TYltpir*  £2  +  2c. 

K:  Iup+(?i}!-lffi)ffr  -±\C  +- 


In  6-1,  Hj  and  Ha-  arc  the  forward  and  alt  propeller  angular  momenta, 
giving  rise  to  the  gyroscopic  moments  indicated.  These  moments  arc  zero 
if  the  forward  and  aft  propeller  speeds  are  the  same  and  the  effective  polar 
moments  of  inertia  are  equal.  Note  that  as  a  consequence  of  the  assumption 
that  the  submarine  is  a  body  of  revolution,  the  ?C  -  ^  axes  are  prin¬ 

ciple  axes  and  all  products  of  inertia  arc  zero. 


6.  Z  TRIM  CONDITIONS  -  TWO  PROPELLERS  OPERATING 
Using  the  set  of  equations  which  define  the  control  forces/moments  (Eq.  5-14  ) 
and  letting  the  cyclic  pitch  components  (  and  )  and  the  change  in 
collective  pitch  (  £  )  be  zero,  the  trim  equations  result:  (the  notation  +/q_ 

designates  the  sign  for  the  forward/aft  propellers) 

f>  4-*J J  (6.2) 

*1*  ‘  sm  K -ia i-C* K&'t.+jCi./ oCjj  (6-3) 

io  3  l/'Oxd-H-,  pur  piM>pi;llcr  {6-1) 

Note  that,  in  straight  and  level  flight,  all  other  components  of  control  force/ 
moment  are  zero. 


At  i  -  l  f»  H  -  V  -  l 


In  order  to  establish  the  trim  conditions  for  the  high-speed  case,  expression 
6-2  Is  equated  to  the  hull  drag.  Thus,  for  balanced  operation  of  the  propellers: 

[C<+  l)]  j./2  ptX'X'l  (6-5) 


A  parametric  investigation  of  the  effects  of  the  variables  Vo  ,  *1*  •  hJ  ,  f\ 
on  propulsion  efficiency  and  other  performance  factors  would  ordinarily  be 
made  in  a  comprehensive  design  study.  However,  no  such  study  was  made 
in  the  present  program.  Instead,  the  value  of  N  suggested  by  Electric  Boat, 
Reference  (10),  was  used  and  a  blade  of  reasonable  else  was  chosen  on  the 
basis  of  rough  horsepower  calculations.  If  the  angle  of  attack  corresponding 
to  maximum  lift/drag  is  inserted  into  equation  (6-5)  a  solution  for  Vs  is 
obtained.  Forward  speed  and  horsepower  then  become  direct  functions  of 
propeller  speed  «Q.  .  Thus  with; 


N 

A 


♦V 

3 


16 

3  ft2 

.  10  rad  (max  / £ 
35* 


and  propeller  speed  3  50  rpm,  the  forward  speed  and  horsepower  are: 


u. 


40.  4fVmc  (2*  knots) 
5740  hp  per  propeller. 


Baaed  on  these  trim  conditions  the  parameters  needed  in  the  computation  of 
the  propeller  derivatives  of  6,4  below  are: 
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F*0  j  */+•  66, l00^ 

Mxo  'V  601,000 #  ft 

7VT  =  y«  2.20 

-  -A  1.75 

"fez.  •=  y+  .211 

The  above  set  of  Constanta  were  used  in  all  of  the  high-speed  stability  and 
control  work  to  be  described  in  Section  VII. 

6.  3  TRIM  CONDITIONS  -  ONE  PROPELLER  OPERATING 
Under  certain  operational  conditions  it  may  be  necessary  or  desirable  to 
operate  the  TPS  with  one  propeller  running  and  the  other  one  locked.  The 
results  of  a  brief  investigation  of  the  required  trim  settings  for  this  mode 
of  operation  are  described  below. 

The  following  principal  assumptions  were  made: 

(1)  The  forward  propeller  is  rotating  in  a  trimmed  condition  at 

some  blade  flight  path  angle  ^  and  angle-of-attack 
corresponding  to  trim  pitch  and  ‘irc 

zero).  The  submarine  is  in  straight  and  level  flight. 

(2)  The  aft  propeller  is  locked  with  the  blades  set  at  some  angle 

of  attack  ^  being  ' 

(3)  The  lift  curve  slopes  for  the  forward  and  aft  propellers  are 
the  same  (i,  e,  ,  the  cascade  effects  are  the  same). 

The  conditions  to  be  met  are  (a)  net  thrust  =  total  drag  and  (b)  the  sum  of  the 
X  -moments  =  zero.  Using  equations  6-2  and  6-3  these  constraints  can  be 
expressed  as: 


A<  1  - 1 6  1 !  -  V  - 1 


/2p  V  "AM  Cos  if  -Sm  5rf  ( Q.  4-  f,  Cl/  <<,*•)] 

-  UpA-XLM ;+-'/*. pc  AM  (C«u  Cl/oC’) 


in  which  the  hull  drag  derivative  is  aaaumed  to  be  valid  at  the  total  speed 


.  and 

-  Vipi^ANK^Xf  SmVf  +  CoS  Vf  [C<a 

+  'Af  K  ANR.CU<^  -  o  <4'7> 


Expression  6-7  can  be  solved  for  the  aft  angle -of-attack: 


<*.  ■  (jf~  r  0/  S«  Yf  +  Co>  (q,  -t-f,  CL/ J/j 

c 

which  is  seen  to  be  a  function  of  ^  and  0^  only  (  a  CoS'SC  ^  ) 

and  not  directly  dependent  on  the  speed  parameters  (  or  tbo  .).  For  any 

particular  set  of  values  )fr  and  oL  ,  o/«.  *s  found  from  6-8. 


If  the  expression  for  o^a  given  by  6-8  is  substituted  into  the  thrust-drag 
equation  6-6  a  quartic  equation  in  the  variable  0^  cast  be  written,  in  which 
the  coefficients  are  functions  of  and  the  physical  constants  of  the  system. 
This  equation,  6-9,  will  contain  within  it  the  thrust-drag  constraint  as  well 
as  the  moment  balance  constraint.  The  expression  is: 


7?  +  AX4  -1-8  -h£  -o  (6-9) 


where 

and 


a  - 

^  -  c.  -Ka  Vp 

e>  =  kf  (i+^n  zh  aa 

c  -  To.v\  ^  ( zc^0f  -  $i<\  Vc-) 

J>  *  ft  f  -fo.*\  iff ) *•  9-  c<o  /.LtllilL)  Me 

(am  Um iff  'HP 
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(Equation  (6-9)  was  solved  (or  various  values  of  between  5*  and  iO*,  I'hc 
values  of  3^  which  are  physically  realizable  and  which  cur  respond  to  them 
^  are  plotted  in  Figure  6-1.  Also  shown  in  the  figure  are  the  corresponding 
settings  of  the  aft  propeller,  found  from  equation  (6-8), 

Having  determined  the  permissible  set  of  jjp  and  ^  which  will  satisfy  the 
thrust-drag  equality,  and  the  corresponding  needed  to  achieve  roll-moment 
balance,  it  is  necessary  to  next  consider  forward  propeller  power.  The 
forward  propeller  power  is: 


P.f  --  M/.„f  -a  -  P„  (  po 


*  O 


(6-10) 


Using  equation  6-3  and  on  nondimensionalizing  (6-10)  by  dividing  by  U,* 

one  obtains  a  nondimensional  power  equation: 

<?•  R  =  hk  t-Cos  (i.„, 

In  Figure  6-2  (p  is  plotted  versus  ^  ,  the  forward  propeller  (blade)  flight- 
path-angle.  A  minimum  occurs  around  23*,  corresponding  to  /  f  of  about 
4.  8°  and  ^  a^out  I’R  7*. 


In  order  to  observe  the  trend  in  the  variation  among  the  variables  l/f  (forward 

o 

speed),  JX.  (propeller  angular  "docity)  and  horsepower,  two  sets  of  calcula¬ 
tions  were  made:  one  at  variable  and  max  /L  (50  rpm),  and  one  it 
constant  and  variable  .fL  .  The  forward  speed  is  common  to  both  sets 

of  calculations.  The  results  (approximate)  are: 


*  Extraneous  values  arise  from  the  squaring  operation  involved  in  eliminating 
and  neglect  of  the  absolute  value  limitation  on  the  induced  drag  for  the 
angle  range  considered. 
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j  Constant  (max) A  ,  50  rpm 

|  Constant  fyp 

li  ,u/L. 

Sf 

Pup 

Vfi., 

if 

rpm 

10.  16 

10* 

478 

.021 

23* 

202  .050 

20.8 

15.5 

15* 

882 

.018 

23* 

720  .022 

31.7 

24.  5 

23* 

2820 

.0087 

23* 

2820  .0087 

50 

27 

25* 

3780 

.0072 

requires  rpm  ^>50 

ser 


Figure  6-2  show*  that  in  order  to  operate  at  maximum minimum 
at  all  speeds,  a  in  the  vicinity  o £  22  or  23  degrees  should  be  used,  and 
this  is  borne  out  by  the  tabulated  data.  Since  =  23*  corresponds  to 
oC^  =  4.8*  and,  since  the  ^  corresponding  to  (max  ^/£)  )*  is  clo 

to  25*,  the  question  arose  as  to  whether  computational  inaccuracies  may  have 
crept  into  Fig.  6-2  to  mask  a  minimum  closer  to  25*.  However,  a  check  of 
the  computations  used  in  this  figure  indicates  otherwise.  Apparently  then, 
optimum  operation  does  not  occur  at  (max  u//[)  ). 

Aside  from  any  question  as  to  what  are  the  exact  optimum  operating  condi¬ 
tions,  it  is  obvious  from  the  above  computations,  that  operation  with  one 

i 

propeller  locked  is  practical  and  that  speeds  of  about  2?  ft/ sec  (16  knots)  are 
possible. 

6.  4  NUMERICAL  DATA 

All  of  the  numerical  data  used  in  connection  with  the  stability  and  control 
studies  are  collected  in  this  section. 

6.  4,  1  Physical  Data 

Physical  data  for  the  hull  and  the  propellers  as  well  as  trim  parameters  for 
balanced  two -propeller  operation  arc  given  in  Table  6-1  for  the  high-speed 
case. 
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6,4.2  Hull  and  Propeller  Hydrodynamics 

The  non-dimensional  stability  derivatives  for  the  hull  (Xh»Y^  ,  etc.)  and 

the  dimensional  derivatives  for  the  propellers  ( * Vv  »  etc.)  are  not 

listed  separately.  Instead,  these  derivatives  are  listed  together,  in  Table 
6-2  where  their  combined  effect  is  accounted  for  by  the  "total"  stability 
derivatives  (Xv*.  t'fv  *  etc* )  noted  earlier  in  Section  V.  In  Table  6-2  the 
relative  effect  of  each  contribution,  hull  and  propeller,  can  be  noted  and 
their  sum  compared  with  the  corresponding  derivatives  of  a  conventional 
submarine.  The  particular  conventional  submarine  for  which  data  were 
available  at  the  time  of  writing  this  report,  is  the  Albacore,  although 
eventually  it  is  hoped  that  comparison  with  a  submarine  of  the  Thresher 
class  will  be  possible. 

Most  of  the  values  for  the  hull  and  propeller  derivatives  given  in  Table  6-2 
can  be  passed  over  without  comment.  However,  there  are  a  few  of  these 
derivatives  which  are  worthy  of  a  few  remarks. 

(1)  Although  at  high  speed  the  contribution  to  ^(g,(rate-of-change- 
of -drag  with  change  in  forward  velocity),  due  to  the  propellers 
is  much  larger  than  that  due  to  the  hull,  the  relative  contribu¬ 
tion  to  total  drag  at  40.4^/  sec  is:  hull  drag^  132,  000#, 
propeller  drag  *=  8100#.  The  propeller  drag,  for  two  propellers 
operating,  can  be  computed  from  the  negative  component  of 

the  trim  thrust  given  in  equation  (6-2). 

(2)  The  propellers  contribute  a  significant  amount  of  damping  in 
heave  (£w  ),  over  that  due  to  the  bare  hull,  but  not  enough  to 
match  the  effect  of  the  empennage  in  a  conventional  submarine. 

(3)  Although  the  propellers  contribute  significantly  to  damping  in 

pitch  and  damping  in  yaw  (  )  over  the  bare  hull,  the 

contribution  is  not  as  effective  as  the  empennage  in  a  conven¬ 
tional  submarine. 
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TABLE  6-2  -  Total  Hull  St  Propellar  Darivativo* 
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6. 4.  3  Control  Coefficient! 

The  numerical  values  for  tho  primary  control  coefficient*,  along  with  the 
corresponding  designated  control  inputs  are  given  in  Table  6-3,  By  way  of 
comparison  with  a  conventional  submarine,  the  Albacore  ^  -force  coefficient 
)  18  about  three  times  as  large,  and  the  pitching  moment  coefficient 
(  JVJj  )  about  four  times  as  large  as  tho  values  in  Table  6-3. 

The  control  cross-coupling  coefficients  are  listed  in  Table  6-4.  The  offect 
of  these  cross-coupling  coefficients  on  control  system  performance  will  be 
discussed  in  Section  VII. 
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TABLE  6-3  -  Primary  Control  Coefficients 
TABLE  6-4  -  Control  Cbupling  Coefficients 
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STABILITY  AND  CONTROL  STUDIES 


7.  1  GENERAL  REMARKS 

The  stability  and  control  properties  of  the  TPS  are  an  essential  element  of  the 
total  number  of  factors  that  bear  on  the  engineering  feasibility  of  the  TPS 
concept.  In  the  most  general  sense,  the  concept  under  investigation  postulates 
a  control  configuration  that  has  greater  control  effectiveness  at  low  (or  aero) 
speed,  plus  greater  flexibility  in  producing  control  forces  and  moments  at  all 
speeds,  than  a  conventional  submarine.  This  gain  in  control  capability  is 
also  accompanied  by  a  loss  in  dynamic  stability  (including  a  loss  in  control 
effectiveness  at  high  speeds)  as  a  result  of  the  elimination  of  the  stabiliaing 
surfaces  (fixed  planes  and  fins)  found  on  the  conventional  submarine.  This 
means  that  the  TPS  will  not  reach  any  attitude,  depth,  or  heading  trim  when 
proceeding  forward  at  high  speeds  unless  the  hydrodynamic  forces  and 
moments  that  cause  the  submarine  to  diverge  from  the  intended  path  are  con¬ 
tinually  corrected  or  counter-balanced  by  control -force  and  moment  inputs. 

This  continual  correction  task  can  be  done  either  manually  by  a  human  con¬ 
troller  or  by  an  automatic  control  system.  If  the  stabilization  task  is  under¬ 
taken  by  the  human  operator,  in  addition  to  his  maneuver-control  task,  the 
magnitude  and  difficulty  of  this  task  enters  into  the  assessment  of  the  handling 


IWife 


qualities  of  the  submarine.  Although  It  may  be  possible  to  stabilize  and  control 
the  TPS  by  manual  means,  the  oxtent  to  which  a  human  operator(s)  could  do 
this  in  acceptable  fashion  Is  open  to  question.  A  glance  at  the  response  of  the 
unmodified  (i.  e. ,  operated  open  loop)  TPS  (see  Figure  (7-1))  shows  that  when 
t(a  -a  40  ft/sec.,  the  pitch  angle,  &  ,  diverges  about  15  degrees  in  about 
15  seconds.  Accordingly,  a  basic  premise  in  this  stability  and  control  investi¬ 
gation  Is  that  automatic-stabilization  methods  must  be  employed. 

The  need  to  employ  control  forces  and  moments  for  stabilization  purposes 
immediately  poses  the  question  as  to  what  extent  the  stabilisation  function, 
whether  done  manually  or  automatically,  will  detract  from  the  maneuverability 
of  the  TPS.  If  part  of  the  available  control  power  must  be  used  to  offset  the 
existing  destabilizing  hydrodynamic  forces,  then  the  total  available  control 
power  can  not  be  used  for  maneuvering  purposes.  In  essence,  this  is  the  cen¬ 
tral  issue  of  feasibility  from  the  stability  and  control  point  of  view,  namely, 
to  what  extent  does  maneuverability  suffer  from  the  requirement  to  stabilize 
with  the  control  function.  Of  course,  there  are  secondary  stability  and  control 
questions,  e.g.,  to  what  extent  are  problems  created  by  the  pitch-yaw  and  yaw- 
pitch  couplings  that  are  present  in  the  cyclic-pitch  control  process. 

The  major  purpose  of  this  initial  feasibility  investigation  is  to  resolve  this 
central  issue  of  stabilization  versus  maneuverability  and  to  consider,  where 
possible,  all  those  characteristics  of  tandem-propeller  control  that  have  a 
bearing  on  the  feasibility  of  designing  an  actual  control  system.  The  overall 
task  remains  to  be  completed  within  the  time  period  of  the  existing  contract. 

To  date,  however,  it  has  been  possible  to  obtain  a  preliminary  reading  on 
whether  the  TPS  has  sufficient  control  power  to  provide  for  both  stabilization 
and  maneuvering  at  high  speeds.  A  comparison  of  the  high-speed  maneuver¬ 
ability  of  the  TPS  with  that  possessed  by  conventional  submarines  remains  to 
be  accomplished,  however. 
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Those  investigations  of  the  dynamic  behavior  of  the  TPS,  completed  to  date, 
are  described  In  the  following  sections.  The  (1)  use  of  linearized  equations 
and  (2)  the  elimination  of  second-order  forces  and  moments  (arising  from 
control  couplings)  facilitated  a  limited  amount  of  theoretical  analysis.  How¬ 
ever,  the  bulk  of  the  stability  and  control  Investigations,  described  herein, 
was  accomplished  with  the  aid  of  analog  computer  simulation.  Representative 
machine  runs  are  included. 


7.  2  THE  LINEARIZED  EQUATIONS  OF  MOTION 

Section  VI  presented  the  dynamic  equations  of  motion  (6-1)  that  describe  the 
six-degree-of-freedom  responso  of  the  TPS  to  control  forces  and  moment*. 
Equations  (6-1)  arc  nonlinear  for  large  motions,  since  they  contain  products 
of  dependent  variables  as  well  as  the  trigonometric  functions  introduced  by 
the  Euler  angle  transformations.  If  small  perturbations  from  trim  are 
assumed,  it  is  legitimate  to  (1)  make  the  small-angle  approximation  that  sin  A  = 
tan  A  s  A  and  cos  A  =  1  and  (2)  neglect  the  products  of  perturbation  variables. 
Equations  (6-1)  then  become  (for  symmetrical  propeller  operation); 
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Note  that  the  forward-velocity  disturbance  variable,  Zt  ,  does  not  cntor  Into 
the  last  five  equations  of  (7-1).  The  latter,  therefore,  define  the  response  of  a 
five-degree-of -freedom  system  in  which  the  forward  velocity,  Ll0  ,  remains 
constant. 

If  we  introduce  the  transformation 

£  -  -  tfosin  £  +  fVcos  Q 

where  J  is  the  vertical  velocity  of  the  submarine  relative  to  an  earth  coordinate 
system  and  then  make  use  of  the  small-angle  assumptions,  there*  is  obtained 

\  -  -  UoO  +  M?  (7-2) 

Substituting  Equations  (5-0)  and  Equation  (7-2)  into  the  last  five  equations  of 
(7-1)  yields: 


r* 
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The  control  forces  and  moments  appear  on  the  right  hand  side  of  Equations  (7-3), 
and  have  been  defined  earlier  by  Equation  (5-15),  If  the  nonlinear,  control¬ 
coupling  terms  are  not  significant  in  the  operating  regime,  Equations  (7-3)  be¬ 
come  a  set  of  linear  differential  equations  and  the  dynamic  behavior  of  the  TPS 
can  be  investigated  using  methods  of  operational  calculus.  Note  that  examination 
of  the  depth-changing  behavior  has  been  facilitated  by  transforming  the  body-axis 
variable,  ,  to  the  earth- refo renced  variable,  • 
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7.  3  ANALYTICAL  STUDIES 

A  brief  analytical  Investigation  of  the  simplified  pitch/ yaw  dynamics  of  the 
TPS  was  undertaken  aa  a  first  step  towards  accomplishing  the  stability  and 
control  work  outlined  above.  Recognizing  that  this  v/aa  to  be  followed  by  a 
more  complete  examination  of  submarine  dynamics  via  computer  simulation 
the  objectives  of  the  study  were  limited  to  (a)  obtaining  an  understanding  of 
the  basic  uncontrolled  behavior  of  the  TPS,  and  (b)  outlining  the  control 
system  requirements  for  pitch,  depth  and  course  control  of  the  vehicle. 

The  uncontrolled  pitch-plane  dynamics  of  the  submarine  can  be  described 
approximately  by  the  and  2,  equations  of  the  set  (7-3).  Focusing 
attention  on  the  response  to  applied  pitching  moment  and  describing  this 
pitching  moment  simply  by  the  following  transfer  function  results: 

(  5  1*  the  Laplace  operator) 

=  _  Ms  _ 

In  a  similar  manner,  the  simplified  yaw  response  to  a  yaw  moment  KI^Sn/ 
can  be  obtained  from  the  Y*  and  N  equations  of  the  set  (7-3).  It  is: 


^  _ _ ~  (»!*>- Yv')  _ _  (7-S) 

If  one  now  substitutes  the  numerical  values  of  the  physical  constants  and 
stability  derivatives  indicated,  each  of  the  resulting  characteristic  deter¬ 
minants  (i.  e.  the  denominators)  of  (  7-4)  and  (7-5)  will  prove  to  possess  a 
negative  coefficient.  Thus,  the  uncontrolled  submarine  is  unstable  in  both 
pitch  and  yaw.  Mathematically,  this  instability  arises  primarily  from  the 
magnitudes  and  signs  of  the  angle-of-attack/  sideslip  derivatives  fAv/  and  A/v 

*  The  word  "uncontrolled"  is  used  to  describe  the  dynamics  of  tl  e  TPS 
without  stability  augmentation  through  automatic  control. 
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but  is  also  due  to  other  dynamic  effect*.  Physically,  It  arise*  from  the 
fact  that  if  $  or  ^  are  disturbed,  a  hydrodynamic  upsetting  moment  rather 
than  a  righting  moment,  is  developed. 

Accordingly,  the  functions  of  any  automatic  control  system  designed  for  the 
TPS  should  be  (a)  to  stabilize  the  inherently  unstable  vehicle,  and  (b)  to 
"shape"  the  dynamic  behavior  of  the  vehicle  to  desired  specifications.  Taken 
together  these  two  functions  constitute  what  is  often  called  "stability  aug- 
mentation'.' 

Based  on  approximations  to  the  higher  order  transfer  function  represented 
by  (7-4)  the  pitch-plane  analytical  investigation  was  carried  to  the  point  of 
determining  crude  ranges  for  stabilization  gains  for  the  pitch  axis.  Since 
this  work  formed  the  basis  for,  and  was  corroborated  by,  an  analog  computer 
study  (see  below)  it  is  not  described  here. 

In  the  sections  to  follow,  stability  augmentation  techniques  which  are  applicable 
to  the  mv  . -degree-of-freedom  control  problem  in  general  are  described 
briefly.  Attention  is  then  given  to  the  pitch-plane  dynamics  (including  an 
analog  computer  investigation).  Finally,  a  few  remarks  are  made  relative 
to  the  yaw-plar.e  dynamics. 


7.4  STABILITY  AUGMENTATION 

Two  stability  augmentation  techniques  were  investigated.  In  the  first,  called 
direct-axis-stabilization  (and  abbreviated  to  D.  A.  S.  for  convenience),  feed¬ 
back  control  signals  are  selected  which  are  functions  of  those  motion  vari¬ 
ables  which  can  be  identified  with  the  instability  of  the  vehicle  or  with  the 
particular  response  which  needs  "shaping"  or  improvement.  These  signals 
are  then  used  to  generate  control  forces  andmoments  in  the  same  axis.  For 
example,  the  term *  which  appears  in  the  pitching  moment 
expression  (Eq.  7-3),  takes  on  a  negative  value  when  .  For 

the  TPS  this  speed  is  about  1.7  ft/  sec.  Above  this  speed  the  vehicle 
exhibits  tho  instability  previously  alluded  to.  In  order  to  stabilize  the  pitch 
axis  let  the  blade  angle  variable  which  produces  pitching  moment,  ,  be 

defined  by: 
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Srv)  =  S/HC  ~  f<A  &  -  K*  9 


(7-6) 


where  i/HC  l*  ^e  commanded  value  of  (1.  o.  the  control -loop  input), 

Sm  if  the  actual  blade  deQection(s)  and  /^|  and  are  pitch-rate 
and  pitch-angle  £u«u'ouck  respectively.  By  proper  choice  o £  the  gains  h(£ 
and  one  can  modify  the  coefficients  in  the  transfer  function  (7-4),  which 
would  then  be  #•  fi).  The  inherent  instability  can  be  removed  and,  within 
limits,  the  dynamic  response  can  be  shaped.  This  effect  can  be  seen  in 
another  way  by  substituting  (7-6)  into  the  A 1  equation  of  the  set  (7-3). 

-AJ7  &  A/* Iff  •  ar  /flj  Svrtc  (7*7) 


in  which  /He  has  been  simplified  to  include  only  the  primary  control  term 
Although  (7-7)  shows  how  the  rate  and  position  gains  and  ^ 
augment  the  damping  and  stiffness  in  the  Af  equation  of  the  uncontrolled 
submarine,  it  does  not,  of  course,  show  the  complete  effect  of  these  feed¬ 
backs  on  the  characteristic  determinant  of  (7-4).  Note  also  that  direct- 
axis  stabilization  does  not  effect  in  any  way  the  hydrodynamic  coupling  terms 
Ai,y  and  in  (7-7).  Coupling  terms  such  as  these  produce  higher- 

ordered  characteristic  determinants,  illustrated  by  (7-4). 

In  order  to  determine  what  benefits,  if  any,  might  accrue  in  reducing  the 
order  of  these  determinants,  a  second  stability  augmentation  technique  was 
investigated  briefly.  In  this  technique  feedback  signals  which  are  generated 
by  motion  variables  in  one  axis  are  used  to  develop  control  forces/ moments 
in  another  axis.  These  feedbacks  can  be  adjusted  to  effectively  cancel  the 
hydrodynamic  coupling  terms  and  for  this  reason  the  method  is  called 
decoupling  stabilization,  A  control  equation  which  illustrates  this  is: 
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which  is  like  (7 -6}  except  that  the  "decoupling  feedbacks"^'  J 
have  been  added.  If  one  now  substitutes  (7-8)  into  (7-3),  as  was  done 
previously  with  (7-6),  it  can  be  seen  that  by  proper  choice  of  and 

the  hydrodynamic  coupling  terms  can  be  made  to  vanish.  It  will  be  shown 
later  that  in  the  case  of  the  pitch  axis,  illustrated  ir.  (7-8),  the  direct-axis 
terms  ^  and  must  be  retained  for  reasons  of  stability.  Hence, 

in  some  of  the  analog  computer  runs  to  follow  this  technique  is  called 
decoupling  -plus-direct-axis  stabilization. 


7.5  PITCH -PLANE  DYNAMICS 

The  analytical  studies,  described  briefly  above,  were  followed  by  an  analog 
computer  simulation  to  further  investigate  the  high-speed  pitch-plane  dyna¬ 
mics  of  the  TPS.  Motion  was  arbitrarily  confined  to  the  plane.  The 
equations  of  motion,  taken  from  (7-3)  are: 


^  ~  /* l (Vfj -h)ko  -  .7  Q  —  2*/ 6 

IiH  9  -  /Wf  6  -  r  /«»•«. 

t 

x,  it  -  Xu. 
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The  and t\  equations  were  mechanized  in  order  to  determine  how  the  non¬ 
linear  control-coupling  terms  contained  in  Xcand  Kc  would  effect  velocity 
and  roll  angle,  A  roll  stabilization  loop  was  used  in  ail  of  the  analog  computer 


*  The  single  and  double  bars  have  no  special  significance. 
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runs.  Thu  roll* rate  and  roll-angle  gains  Kj  and/flj  (see  below)  were  arbi¬ 
trarily  selected  to  yield  an  undamped  natural  frequency  of  about  1.0  rad/  sec 
and  a  damping  ratio  of  0.7  critical  for  the  uncoupled  roll  dynamics. 

Thu  control  equations  are: 
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Both  direct-axis  and  decoupling  gains  are  shown  in  the  equations,  however 
not  all  of  these  were  used  in  every  analog  run.  Note  that  commanded  values 
ate  differentiated  from  actual  values  by  the  subscript  C  . 
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Tho  control  forced  and  momenta,  with  motion  confined  to  the  plane,  are 
given  by: 
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7.5.1  Discussion  of  Analog  Runa  -  High-Speed  Cti«c 

The  reaponae  of  the  basic  submarine  with  all  atability  augmentation  terma 
set  to  zero  and  no  control  coupling  terma  includod  ia  8hown  in  Figure  7-1. 

In  all  of  the  analog  computer  work,  runa  were  made  with  and  without  the 
nonlinear  components  of  (7-11)  and  the  effect  of  these  components  was  found 
to  be  negligible.  These  terma  were  therefore  left  out  of  the  runs  reproduced 
in  thir  report  (except  for  Figure  7-6).  Time  histories  of  pitch  angle,  depth 
change,  roll  angle  and  forward  velocity  are  shown  in  Figure  7-1,  reading 
from  top  to  bottom.  The  instability  of  the  uncontrolled  submarine  is  obvious. 
Notice  that  roll  angle  and  forward  velocity  perturbations  are  too  small  to 
appear  in  the  traces.  In  all  runs,  ^  and  M,  were  monitored  and  found 
to  be  negligibly  small. 

In  Figure  7-2,  feedback  gains  necessary  to  decouple  the  pitch  and  depth 
responses  wore  used.  The  response  to  a  SmC  input  (equal  magnitudes  of 
Xifc  and  <f/&&  )  is  a  divergent  pitch  angle,  but  the  depth  remains  constant. 
The  pitch  terms  in  the  depth  equation  have  been  effectively  cancelled  by  the 
decoupling  feedback.  However,  the  system  is  still  unstable  in  pitch.  In 
this  figure,  and  in  all  of  the  remaining  analog  runs,  the  forward  and  aft 
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propeller  anglen-o£-attack  arc  traced.  These  are  tho  total  change*  In  U.  , 
over  and  above  the  trim  c /  ,  <:,i.U**d  for  .is  a  result  o?  the  particular  feed¬ 
backs  being  used. 

In  Figure  7-3,  direct-axis  feedback  is  added  to  the  system  of  Figure  7-2. 

The  gains  were  selected  to  ahapo  the  pitch  natural  frequency  to  0.  1  rad/ sec, 
with  0,  7  critical  damping.  There  is  no  depth  change  because  ©  • 

Note  that  theo<^  needed  to  achieve  a  steady- state  pitch  angle  of  about  4* 
is  approximately  25*.  This  ar.gic  may  be  beyond  the  actual  stall  angle-of- 
attack  of  the  blades  but  this  depends  upon  blade  design.  In  any  event  blade 
angles-of-aitack  in  excess  of  20  or  25  degrees  were  considered  undesirable 
for  the  purposes  of  this  study. 

In  Figure  7-4,  direct-axi  i  pitch-rate  stabilisation  only  is  used,  (pitch- 
angle  gain,  The  £  trace  of  previous  runs  has  been  replaced  by 

the  -  ^  trace.  A  steady  pitch  angle  of  about  9  degrees  and  a  climb- rate 
of  about  7  ft/eec  are  obtained  after  about  500  sec,  fot  »iep  input  of 

.05  rad.  The  angles -of-attaek,  which  are  at  the  commanded  value  at  f*  C  , 
rapidly  reverse  as  the  vehicle  pitches  up  (the  effect  of  the  feedback) 

and  eventually  settle  out  to  about  .  025  rad. 

The  relatively  poor  response  shown  in  Figure  7-4  can  be  improved  by  adding 
pitch-angle  feedback,  q  .  This  improvement  is  shown  in  Figure  7-5, 
where  it  is  seen  that  the  pitch  and  depth- rate  responses  reach  95%  of  their 
final  values  in  about  SO  seconds.  The  run  of  Figure  7-5  is  repeated  in 
Figure  7-6  with  the  nonlinear  control -coupling  terms  c*f  equation  (7-11) 
included  in  the  simuiuti  »u  in  order  to  show  that  their  effect  is  negligible. 

The  response  of  the  TPS  to  input,  when  direct-axis  stabilization  is 

employed,  .s  shown  m  Figure  7 -7.  ru  this  case,  th0  submarine  does  not 
reach  a  steady- state  condition  where  th-*  hull  anal*  of  attack  is  zero  (as  is 
the  ca»e  in  Figutes  7-5  and  7-0),  unless  Li  Ms 
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FIGURL  7-4  Direct-Axis  St.ibihzal  ion  (Xu  Control  Coupling)  -  Linear 
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Tho  response  times  for  Figures  7*6  and  7*7  are  about  the  name,  but  the 
steady- atato^a  re  higher  and  tho  depth-rate  is  lower  in  Figure  7-7.  In  othc*" 
words,  as  should  be  anticipated,  more  efficient  diving  or  climbing  is 
obtained  with  pitching -moment  input  than  with  £  -force  input. 

Combined  direct-axis  plus  decoupling  stabilization  is  explored  in  Figure  7-8, 
where  a  depth-rate  time  constant  of  14.  3  seconds  and  a  pitch-angle  natural 
frequency  and  damping  ratio  of  .  10  rad/ sec  and  0.70  respectively  are 
employed.  The  decoupling  and  direct-axis  feedback  sensitivities  used 
result  in  steady-state  pitch  angle  and  depth  rate  responses  to  a  slep  input 
given  by  the  following  equations; 
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(7-12) 

-  V 

(7-13) 

In  order  to  cause  the  submarine  to  point  in  the  direction  of  its  total  velo¬ 
city  vector,  both  7mc  and  command  inputs  must  be  used  simultaneously. 

It  can  be  shown  that  having  ^  *  —  U0  Q  (aero  hull  angle-oi-attack)  requires 
that  SjC Note  that  £jc«  •!*  and  that  J z-  4o&  in  Figure  7 -8. 

The  blade  angles-of-attack  are  very  small  and  the  response  times  are  fast 
for  this  case. 


Figure  7-9  indicates  that  responses  similar  to  those  of  Figure  7-8  can  be 
obtained  by  direct-axis  stabilization  alone  and  through  the  use  of  pitch-angle 
command  inputs  alone.  Response  times  of  about  thirty  seconds  result  from 
the  gains  indicated  on  the  figure.  These  response  times  can  be  reduced 
even  further  by  increasing  the  }\ q  feedback  gain.  It  can  be  concluded  from 
Figures  7-8  and  7-9  that  either  stabilization  technique  can  be  used  with 


similar  rosults  in  terms  of  pitch  angle  and  depth  rate  dynamic  and  static 
responses,  and  that  about  the  same  "economy"  in  using  up  the  useful  range 
of  blade  angle -of-attack  is  achieved.  Accordingly,  there  seems  to  be  little 
to  recommend  the  decoupling -stabilization  technique,  in  the  pitch-plane, 
when  climb  or  dive  maneuvering  is  involved,  since  a  £  sensor  would  be 
needed  to  implement  m  the  control  equation  (7-10).  As  has  been 
shown  above,  climb  and  dive  maneuvering  can  be  performed  efficiently 
with  direct-axis  feedback  and  i  inputs  alone.  However,  in  the  case  of 
an  attitude  stabilizing  a.td  depth-keeping  control  system,  or  in  the  case 
of  docking  maneuvers,  there  may  be  some  advantage  to  decoupling  stabiliza¬ 
tion.  This  matter  will  be  explored  further  in  the  work  that  remains  to  be 
done  under  the  present  contract. 

7.6  YAW-PLANE  DYNAMICS 

At  the  time  of  writing  this  report  the  yaw-plane  analog  computer  studies 
were  still  in  orogress.  However,  preliminary  calculations  have  shown 
that  the  stabilization  problems  in  the  yaw  plane  are  quite  similar  to  those 
in  the  pitch  plane.  It  should  be  possible,  therefore,  to  achieve  stable  and 
acceptable  controlled-submarino  behavior,  but  a  firm  conclusion  to  this 
effect  cannot  be  made  at  this  time.  The  outcome  of  the  analog -computer 
study  will  indicate  whether  or  not  the  strong  yaw- roll  coupling,  due  to  the 
presence  of  the  sail,  is  critical  in  terms  of  the  demands  made  upon  the 
available  control  forces  and  moments.  Thus,  while  it  is  a  relatively  simple 
matter  to  calculate  the  stabilization  feedbacks  needed  to  achieve  specified 
yaw-roll  dynamics,  the  question  is  one  of  whether  these  stabilizing  signals 
will  call  for  more  control  power  than  is  available. 
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VIII 

REMAINING  STUDIES 


The  work  scheduled  in  the  present  program  is  aimed  at  performing  or 
completing  the  following  tasks: 

(1)  Complete  the  low-speed  hydrodynamic  studies, 

(2)  Complete  the  analog  computer  study  of  the  high-speed  yaw- 
plane  dynamics. 

(3)  Perform  analog  computer  studies  of  the  low-speed  controlled- 
submarine  dynamics. 

(4)  Draw  preliminary  conclusions  as  to  the  feasibility  of  the 
tandem  propeller  concept  based  upon  the  hydrodynamic  studios, 
the  stability  and  control  studies,  and  judgments  (yet  to  be  made) 
of  the  estimated  handling  qualities  of  a  Tandem  Propeller 
Submarine. 


HH 


All  - 1 6  H  •  V  - 1 


,nr,  ,«<i  '«■»  fi; •'  HjR»  ftf  A 


l 


IX 

REFERENCES 


1.  "Naval  Architectural  Aspects  of  Submarine  Design"  by  S.  S.  Arcntzen 
and  Philip  Mandel.  Trans.  Vol.  68,  19t>0  of  the  Society  of  Naval  Archi¬ 
tectural  k  Marine  Engineers  (U). 

2.  "Qualitative  Discussion  of  the  Hydrodynamic  Properties  of  Tandem  Propeller 

Systems  With  Cyclic  Pitch  Control"  Netherland  Ship  Model  Basin,  February 
1962  <U). 

3.  "Force  and  Pressure  Distribution  Measurements  on  Eight  Fuselages"  by 
G.  Lange.  NACA  TM  #1194,  October  1941. 

4.  "Model  Investigation  of  the  Stability  and  Control  Characteristics  of  the 

♦ 

Preliminary  Design  of  the  SSGN  (FBM)  by  D.  B.  Young,  David  Tayl<  r 
Model  Basin,  Report  #C-928,  July  1958. 

5.  "A  Method  for  Calculation  of  Hydrodynamic  Lift  for  Submerged  and  Planing 
Rectangular  Lifting  Surfaces"  by  K.  L.  Wadlin  and  K.  W,  Christopher. 

NASA  TU -It -l 4.  1959. 


H‘)  AC  i  -  I  n  i  I  -  V  - 1  f 


AG  - 1  o  11  -  V  - 1 


6.  "The  Stability  Derivatives  of  an  Airscrew"  by  H.  Glauert,  Koyai  Air¬ 
craft  Establishment  K  k  M  #644,  1919  (U). 

7.  "Aerodynamics  of  the  Airplane"  (book)  C.  B.  Millikan,  John  Wiley  &  Son 
1941. 

8.  "Airfoil  and  Airscrew  Theory"  (reprint)  by  H.  Glauert,  Cambridge 
University  Press  1943. 

9.  "Large  Hub  to  Diameter  Ratio  Propeller  With  Programmed  Blade  Control" 
Netherlands  Ship  Model  Basin,  January  1962. 

10.  "Feasibility  Study,  Tandem  Propeller  Propulsion  System  for  a  Submarine, 
Electrical  Machinery  Specification,  Revision  C"  U413-61-192,  December 
19,  1961.  (U). 


90 


AG  - 1 6  S  \  -  V  - 1 


APPENDIX 

THE  LINEARIZED,  LOW-SPEED,  HYDRODYNAMIC  FORCES 
AND  MOMENTS  ARISING  FROM  THE  PROPELLERS  OF  THE  TPS 


The  force  and  moment  equation*  resulting  from  the  analysis  described  in 
Section  5.  2.  2.  2  are  presented  in  this  appendix.  The  equations,  as  shown, 
are  valid  for  the  forward  propeller  only  when  this  propeller  is  producing  a 
positive  thrust.  Modifications  in  the  sign  of  various  terms  are  required  to 
make  the  equations  valid  for  the  aft  counterclockwise  rotating  propeller  and 
to  further  handle  those  situations  where  the  submarine  axial  motion,  CL  ,  is 
opposite  to  the  direction  of  the  propeller-induced  inflow  velocity,  .  It 
should  be  noted  that  the  equations  are  based  on  body  axis,  and  that  all  terms 
involving  (  a  )  and  (  )  have  been  multiplied  by  2  to  account 

for  the  cross-flow  velocity  distribution  around  the  hull. 
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